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A unitary multi-channel approach, first developed by the Carnegie-Mellon-Berkeley group, is 
applied to extract the pole positions, masses, and partial decay widths of nucleon resonances from 
the partial wave amplitudes for the transitions from nN to eight possible final baryon-meson states. 
Results of single energy analyses of the VPI group using the most current database are used in 
this analysis. A proper treatment of threshold effects and channel coupling within the unitarity 
constraint is shown to be crucial in extracting resonant parameters, especially for the resonance 
states, such as Sii(1535), which have decay thresholds very close to the resonance pole position. 
The extracted meisses and partial decay widths of baryon resonances up to about 2 GeV mass are 
listed and compared with the results from previous analyses. In many cases, the new results agree 
with previous analyses. However, some significant differences, in particular for the resonances that 
are weakly excited in nN reactions, are found. 
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I. INTRODUCTION 

The interest in the study of baryon resonances began many years ago and led to the important discovery of SU(3) 
symmetry. Although many states were discovered, the quality and scope of the data limited the analyses. Interest 
has grown significantly in the last few years because of the prospects for new data of high quality at facilities such 
as the Thomas Jefferson National Accelerator Facility, the Bonn synchrotron, and Brookhaven National Laboratory. 
These are excited states of the nucleon with excitation energy of about 0.3 - 1 GeV. Since most of these states were 
discovered in amplitude analyses of ttN —> nN scattering data, they are labeled by the approximate mass and the ttN 
quantum numbers: the relative orbital angular momentum L, the total isospin T, and the total angular momentum 
J. For example, the I?i3(1520) resonance has a mass of about 1520 MeV, isospin T=l/2, total angular momentum 
J=3/2, and decays into a L=2 ttN state. Information about the various baryon resonances is evaluated and tabulated 
in biennial publications of the Particle Data Group |^]. In their 1998 publication, they list about 20 N* (T=l/2) 
and 20 A (T=3/2) states. Some states are well established in the data and most analyses agree on their properties. 
On the other hand, some very large discrepancies exist between different analyses. For example, the extracted full 
widths for the Sii(1535) state (A 4* state of PDG 0) are 66 MeV @, 120±20 MeV §, 151±27 MeV [|, 151-198 
MeV 1^, and 270±50 MeV ||]. This state has a number of striking properties - unusual strong decay width to r]N 0] 
and an unusually flat transition form factor 0] - that have made it a very interesting state to understand. Although 
there are a number of strongly excited states (rated 4* or 3* by PDG), there are many weakly excited states whose 
properties are poorly determined with the existing data. The differences between various amplitude analyses originate 
for various reasons, including handling of data, method of parameterizing nonresonant background, and handling of 
the many N* decay channels. We emphasize the fact that the baryon resonances always decay into a baryon and one or 
more of various mesons; the existing analyses differ from each other significantly in the methods used to describe this 
intrinsically multi-channel problem. Thus, model dependence in extraction of resonance properties makes it difficult 
to test predictions of theoretical models with existing data. With the new experimental facilities, the situation will 
soon be greatly improved when more exclusive data for different final states, such as rjN, ttA, and pN, become 
available. One of the major goals of the new experiments is to obtain data for a large number of reactions. This 
greatly increases the probability of seeing new resonances, e.g. those that couple weakly to the ttN channel, but 
puts significant demands on models to interpret the data consistently. The objective of this work is to revive the 
multi-channel analysis of partial wave amplitudes into resonance parameters of the Carnegie-Mellon Berkeley (CMB) 
group. Although this is meant to be a significant step toward a full analysis of all contributing reactions, the equally 
important problem of extracting partial wave amplitudes from the observables is not considered in this paper. 

The baryon resonances (called N* states from now on in this paper) extracted from amplitude analyses are thought 
to be predominantly composed of 3 valence quarks because of the SU(3) symmetry seen in the spectrum of hadrons 
of low total angular momentum. This notion has been the basis for developing various quark models, ranging from 
the well-studied Constituent Quark Models ||-|ll|, Chiral Bag Models [|l2|,|^, NJL models @, Sohton models [^, 
to the most recent Chiral Constituent Quark Models jl^l. All of these hadron models are motivated by QCD and 
are constructed phenomenologically to describe the hadron spectrum in terms of suitable parameterizations of quark 
confinement mechanisms and quark-quark interactions. With some additional assumptions on the decay mechanisms, 
these models can use the constructed wave functions to predict the decay widths of the N* states. For example, the 
one-body form of quark currents is used ||l7|, [l0| i n most of the constituent quark models in predicting the N* — > ^N 
decay width. A ^Pq model is assumed in Ref. |18| for creation of qq pairs in the calculation of the decay of an N* into 
nucleons and mesons. The comparison of these predictions with the decay width data, such as that listed by PDG, 
is clearly a more detailed test of hadron models. In recent years, lattice QCD calculations have been extended to 
predict masses of the low lying baryons. iQ. Although the lattice QCD calculations are far superior to the empirical 
quark models, they are much more difficult to apply to the decay widths. In any case, precise data of resonance 
parameters such as masses and partial decay widths for various final states are needed to distinguish QCD-inspired 
hadron models and, ultimately, understand the non-perturbative aspects of QCD. 

The determination of resonance parameters is almost always a two-step process. First, phase shift analyses (elastic 
data) and isobar analyses (inelastic data) are used to separate the cross section and polarization observable data into 
partial wave amplitudes, often in the form of T matrices. With a large data set, this determination has less model 
dependence than the extraction of resonance properties in the second step. If nonresonant effects are small, resonances 
will show up as counterclockwise rotations in the complex energy plane (Argand diagram) as the energy dependence 
of these T matrices is plotted. Partial wave amplitudes have been determined by various groups for elastic scattering 



and much less often for inelastic scattering. Older elastic scattering analyses of Carnegie-Mellon Berkeley (CMB) EG] 
and Karlsriihe Helsinki (KH) l^^^fj] stressed the importance of theoretical constraints such as dispersion relations to 
ensure the uniqueness of the fit. The data situation has improved significantly since then; the more recent work of 
the Virginia Polytechnic Institute and State University (commonly called VPI) group is the most visible recent effort 
B. They have regularly updated the data base and attempted to cull out older less viable results. They determine 
single-energy (also called energy independent) T matrices for partial waves up to L=6 and satisfy a smaller set of 
dispersion relations than the earlier work. All these analyses are available in the VPI repository 0]. We use the 
VPI single-energy elastic partial wave amplitudes from the 1995 analysis (SM95) in this work. At the same time, we 
recognize the point made by Hohler [^ that this analysis is not as well constrained as the older analyses. We use 
the single-energy solutions rather than the more debatable smoothed solutions. 

The second step is to extract the resonance parameters from the partial wave amplitudes. In a simple picture, 
each transition amplitude in a ttA'^ reaction would be parameterized as the product of the excitation strength of the 
incident channel to a given resonance and the decay strength of the resonance into the allowed final states with a 
resonance propagator for the intermediate state. However, this s-channel resonant mechanism is far from complete 
since the u-channel and t-channel mechanisms, as implied by crossing symmetry or meson-exchange mechanisms, are 
known to be important. Thus the parameterization of the amplitude in each partial wave must contain a resonant 
part and a non-resonant part (called the background term in most of the literature). Furthermore, the threshold 
effects associated with each decay channel(7riV, rjN, jN, ttA, pN, uiN, ttN* (1440) and others) must be treated correctly 
within the multi-channel unitarity condition. Thus, resonance extraction requires a significant calculational effort and 
many articles have presented various ways to proceed in practice. The PDG mostly bases its recommended values 
for baryons on a few works that study the full resonance region (1170 MeV< W <2200 MeV). These include older 
work by the CMB group |p| and the KH group , and more recent work by the Kent State University (KSU) group 
0, and the VPI group ^. A very recent work of Feuster and Mosel [|| fits data for W <1900 MeV. All of these 
efforts use the data of ttN reactions. All maintain unitarity, though the methods employed are quite different. This 
is reasonable since there is more than one way to implement the unitarity condition. However, these analyses differ 
significantly from each other in handling the multichannel character of the ttN reactions. The CMB and KSU groups 
use a formalism that allows for many channels, while the KH and VPI groups focus on the ttN elastic channel. KH 
accounts for all inelasticity in absorption parameters and VPI uses a dummy channel to account for all inelasticity. 
Feuster and Mosel B fit elastic scattering and inelastic scattering cross-section data with asymptotic two-body final 
states directly, and account for all the remaining components of the total inelastic cross section with a dummy 
inelastic channel. For most strongly excited, isolated states, these five analyses tend to agree within the assigned 
errors. However, as mentioned above, significant differences between them exist in many cases. 

This paper will revive the CMB approach [g| and apply it to extract baryon resonance parameters from partial 
wave amplitudes of ttN reactions with a large variety of final states. This approach emphasizes the analytic properties 
of scattering amplitudes in the complex energy-plane that are consistent with the dispersion-relation approach and 
potential scattering theory. Other methods of handling multi-channel unitarity are through the iiT-matrix approxima- 
tion [ppi^jsll and the KSU model [Q. Although the transition amplitude is parameterized in a form similar to these 
models, there are distinct and important differences. In particular, the threshold for each of eight possible channels is 
treated correctly with 2- and 3-body unitarity requirements imposed. Thus, resonances can be found as poles in the 
T-matrix, a feature missing in most A'-matrix models. Obviously, the CMB approach is most suitable for extracting 
resonances that are close to inelastic channel thresholds. Although the KSU approach by Manley and Saleski and 
more recent K matrix models H also account for the structure due to channel openings, their multi-channel pa- 
rameterization is completely different from the CMB model in realizing the unitarity condition. In the CMB model, 
dispersion relations are used to guarantee analyticity in the amplitudes. The KSU parameterization is an extension of 
a A'-matrix formulation. It does not allow the analytic continuation into the complex energy plane, which is required 
to find the resonance pole. Feuster and Mosel present poles from a speed plot analysis, stating that technical problems 
associated with their partial wave decomposition prevent direct determination of poles from their K matrix analysis. 
There is the additional problem that resonances are associated with poles in the T matrix rather than the K matrix. 

In the CMB model, the non-resonant t-channel and u-channel mechanisms are simulated by including the transitions 
to sub-threshold one-particle states with masses below the nN threshold and with the possibility of producing either 
an attractive or a repulsive potential. Existing models treat these effects in a variety of ways. The empirical method 
used here is rather different from the polynomial parameterization often used in previous amplitude analyses. The 
consequence of this difference is very significant in practice since the intrinsic "structure" of the non-resonant term due 
to the channel opening is built in correctly in the CMB approach, but can be easily missed if the "smoothness" of the 
background term is the only criterion. Modern K matrix methods g fit coupling strengths in various diagrams. These 
have the advantage of fitting fewer parameters than the more empirical approaches and have the correct partial wave 
decomposition of the diagrams include. The disadvantages include ignoring all off-shell intermediate state scattering 
and the inability to describe resonances of total angular momentum larger than 3/2. 



The CMB approach was pubHshed in 1979. Although inelastic data were used in the analysis, the elastic ttN 
scattering data were emphasized. In this work, we follow closely their theoretical approach, making only small 
changes in the parameterization of the amplitudes. However, we make significant changes in the data set used. We 
include all TriV -^ rjN data, some of which was unavailable then. (They depended on a separate analysis of backward 
■K~p elastic scattering and other data [p3 to model the eta cusp. This analysis would give different results with the 
present data set.) We also directly fit TriV ^ tittN inelastic data, represented by the amplitudes determined in the 
isobar model fit of Manley, Arndt, Goradia, and Teplitz p^. By including 30 percent more data, the Manley et al. 
inelastic amplitudes are more accurate than those used in the original CMB analysis. Furthermore, we will use the 
most recent VPI energy independent amplitude ^ as the TriV elastic scattering input. This amplitude is significantly 
different from the one used by the CMB group and is also more accurate than the KH amplitude in the Sn channel, as 
discussed in Ref. [^. These differences in the input data make our results for some cases (in particular the Sii(1535) 
state) significantly different from the CMB values listed in Particle Data Table. 

Batinic et al. p7]| have applied the CMB model to perform an analysis with only 3 channels: ttN, r/N, and a dummy 
channel meant to represent the complex set of nnN channels. Their focus is on the dynamics associated with the rjN 
channel. They fit the KH80 [^ energy-independent amplitudes of the nN elastic data and the ttN — > rjN data. The 
rjN -^ rjN amplitudes are the predictions of the model. The Brown et al. data [Esf is the largest body of TriV -^ rjN 
data in the energy region close to threshold, but it is felt to have significant systematic errors [ p7|j29| . Batinic et al. 
give it a weight 5 times smaller than the other data to deemphasize it. We will use the same ttN —> rjN data with 
the same weights in our analysis. This is an important part of the determination of the parameters associated with 
the Sii(1535) resonance. 

The present work is the first step in an ongoing program to develop a model appropriate for the new generation of 
data coming from Jefferson Lab and other labs. Although we anticipate further development of the model, the goal 
of this work is to apply the CMB model in a form very similar to its original implementation. We closely follow the 
procedures of Ref. [|6| , but use it to determine the baryon spectrum and decay branching fractions from the modern 
data. Future publications will address interesting but nontrivial issues such as how to determine the "best" baryon 
spectrum and how to differentiate between resonance poles, bound state or virtual poles, and poles that have changed 
Riemann sheet by moving across an inelastic threshold cut to the sheet most directly reached from the physical data. 

In section II, we present a detailed account of the CMB multi-channel unitary model. To illustrate the main features 
of the model, we discuss in section III the examples of an isolated resonance and the 2 resonance-2 channel situation. 
The full results are presented and discussed in section IV. A summary and outlook are given in section V. A more 
complete discussion of the methods and results is given in Ref. Q . 

II. THE CMB UNITARY MULTI-CHANNEL MODEL 
A. Representation of 3-body final states 

The first step of a multi-channel analysis to determine N* parameters is to extract from the nN — > ttttN data a set 
of partial- wave amplitudes for the transitions from a 'kN state to various quasi two-particle channels in which one of 
the two particles is either a ttN or a tttt resonant state: ttA, pN , 7rA^*(1440) and a very broad {■k-k)j=t=oN channel. 
For convenience, we will use the notation a to label the S-wave, isoscalar tttt state. This procedure was introduced 
by the Carnegie Mellon-Berkeley (CMB) group in the 1970's. A similar procedure was used later by the Virginia 
Polytechnic Institute and State University (VPI) and Kent State University (KSU) groups |^ to obtain partial wave 
amplitudes from all of the vrA^ -^ ttttN data available in 1984. These amplitudes, called VPTKSU amplitudes, will 
be used in this work. 

A more precise analysis should start with the original data of t:N -^ ttttN reactions. This event data has been 
stored by the VPI group. The raw t:N — > nnN data was not received from them in time for the present analysis. 

B. Model details 

To extract resonance parameters from VPTKSU amplitudes, it is necessary to employ a multi-channel formulation 
of the ttN reaction. This is accomplished in the CMB approach by assuming that the transition amplitudes of the 
ttN reaction can be written in the center of mass (cm.) frame as 
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Tab= X! fa{s)^/ Pa{s)laiGij{s)^jb^/ Pb{s)fb{.s) (1) 



where s is the total center of mass energy squared, indices a, b denote the asymptotic channels which can be either a 
stable two-particle state {ttN, rjN, KA, ••) or a quasi two-particle state (ttA, pN ■ •). These asymptotic channels are 
coupled to a set of intermediate states (resonances) denoted by indices i,j. The scattering matrix defined by Eq. dl^) 
is related to the S'-matrix via S" = 1 -f 2iT with S^ S = 1. Hence we have the following unitarity condition: 

M 

Qm{Tat)=J2T:cTcb (2) 

c=l 

The crucial step of the CMB model is to choose a parameterization of various quantities in Eq. (^ such that 
Eq. (Q) is satisfied. Furthermore, the resulting analytic structure of the scattering amplitudes is consistent with the 
well-developed dispersion relations for TriV elastic scattering and multi-channel potential scattering theories. This is 
accomplished by using the following prescription. First, it is assumed that the jth resonance to be found is identified 
with the bare particle, i, with a bare mass squared, sq.j. The strength constant jat and form factor fais) in Eq. (|l|) 
define the decay of the ith resonance into an asymptotic channel a. The form factor is defined by 

/.(.) = (- ^P==] (3) 

where Qi and Q2 are empirical constants defining how quickly a channel, with orbital angular momentum la, opens 
up. Pa is the center of mass momentum for channel a. In this analysis, we set Qi and Q2 equal to the pion mass. 
These are the same values as were used by CMB; they chose this parameterization and these values as a result of a 
study of the best way to model the left-hand cut. 

The right-hand or unitarity cut is set through the channel propagators, 3?7i0a(s) — faPa{s). The phase space factor 
Pa in Eq. (1) is defined by 

Pa = ^= (4 

for a stable two-particle state. For quasi two-particle channels, pa clearly must be defined consistently with the phase 
space factors used in defining the propagation of this resonant two-body state during the collision, as required by the 
unitary condition Eq. (2). In the CMB model, this is achieved by assuming that the only interaction during collisions 
is a vertex interaction which converts intermediate states into asymptotic states. Then the propagator Gij in Eq. dll) 
can be graphically depicted in figure n^ and is defined by the Dyson equation: 

Gij = ^y + Gil^lkGkj (5) 

Eq. (ph is an iterative equation. A sum over repeated indices is implicit. G, Gq, and E all vary with s and are 
N X N matrices, where N is the number of intermediate resonant and nonresonant) states considered. Each bare 
intermediate state has a propagator. Go, defined by: 

Gl{s) = ^M^ (6) 

where e^ — +1 for states that correspond to the resonances that will be fit. To simulate the i— channel and m— channel 
mechanisms, two subthreshold bare states with a mass so,i below the ttN threshold are introduced. These subthreshold 
states will simulate an attractive background potential for e^ = —1, and a repulsive potential for e^ = -1-1. In principle, 
this prescription can simulate any t— and m— channel mechanisms if a sufficiently large number of bare states are 
included. This is well known in potential scattering theory. The actual number of the subthreshold states needed is 
not intrinsically known. We use two subthreshold states and one state at very high energy for every partial wave. In 
the original work, CMB Q found negligible differences between fits using a (smaller) number of subthreshold states 
and fits using actual potentials that simulate the left-hand cut. They are allowed to couple to TriV and rjN asymptotic 
states. 

The self-energy, Y,ij, describes the dressing of bare particles by the coupling with two-particle channels, as depicted 
in figure |[ It therefore must depend on the strengths(7ia's) and form factors(/(i's), and is assumed to take the 
following form 
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5^u(s)=^7cz$c(s)7cj (7) 

c=l 



For the contributions from stable two-particle channels, we have 

Ms) = Ms) (8) 



with 



1 r°° 

Ms) = z ds/;(s')ffc(s';s)/c(s') (9) 

' Sth.c 



TT 



Here Sth.c is the threshold for channel c, and gds) is the propagator of the two-particle channel c. The task now is to 
choose gds'] s) such that the unitarity condition Eq.(2) is satisfied and the desired analytic structure of the amplitude 
can be generated. It can be shown that the (f>c{s) of the CMB model is obtained by assuming that 



is' -,8) = - 

n 



Pcjs') p Pcjs') 



So 



(10) 



where the density of states Pc(s') has been defined in Eq.(4), and P means taking the principle-value part of the 
propagator. Substituting Eq.(lO) into Eq.(9), we then obtain the following dispersion relation for the auxiliary 
function (f>c 

QmiMs))^f!is)pcis) 
^e (Ms)) = ^e (Mso)) + '-^ C W^?0^)ds' ^'^^ 

We see that Eq.(ll) is a subtracted dispersion-relation which has a form similar to what has been established in many 
studies of ttN scattering. In the complex-s plane, one can choose the subtraction point sq such that the resulting 
scattering amplitude has a pole on the left-hand side and a branch cut from sth,-nN = {rriTT + tun)'^ to -|-cxd. For each 
channel, we choose the value at threshold for sq, and we set the subtraction constant, 5Re {4>{sq)), such that 5Re (0c(s)) 
is at threshold. This arbitrary choice has the primary effect of shifting the value of the bare mass squared, Si, for 
each resonance in a partial wave, but does not affect the physical mass of the "dressed" resonance. Here we assume 
that this dynamical assumption is valid for all stable two-particle channels like rjN , and KK. 

For a quasi two-particle channel c, the function $c('5) in Eq.(7) must account for the mass distribution of one of 
the two particles which is itself a resonance state. To be specific, let mi be the mass of the stable particle in the 
channel c and TO2 and 7713 be the masses of the two daughter particles from the decay of the resonant subsystem into 
a channel r. Then the form assumed by the CMB model can be more explicitly written as 

/•(Vs-mi)^ 
$c(s) = / dSr(T[Sr)(i)c{Sr) (12) 



where the mass distribution of the quasi particle was taken to be: 

^ >Sm(c^,(.))/7r 

(Af2-s)2+72CJ„j(^^(5))2 

Here (j)r{s) is again defined by the subtracted dispersion relation Eq.(ll) for the appropriate resonant subsystem. 
Eq.(n3[) has the commonly used Breit-Wigner resonance form. The coupling strength 7^ is related to the width F^ 

of the considered resonance state by Tr = — — jj — ^. In this work r is either a -kN state or a tttt state. From the 
empirical values of the widths for A and p, mass distribution functions for the ttA and pN channels can be fixed. For 
the other quasi two-particle channels, aN and 7riV*(1440), the width F^ is fixed at a standard value in the fit Q. The 
above formalism for the quasi two-body channels provides a unitarity cut along the real s-axis from the three-body 
TTTriV threshold to infinity. Furthermore, it makes the resulting scattering amplitudes satisfy the unitarity condition 
Eq.(2). This is how the three-body unitarity is implemented in CMB model. 

The Dyson equation, Eq.(5), is algebraic and can be solved by inverting a N x N matrix. Schematically, we have 



G^ds) ^ [H-\s)]., (14) 



with the matrix element of H defined by 



Hi]{s) = - ^^%-Sy(s) (15) 

Now all of the ingredients needed for calculating the T-matrix elements of Eq.(l) are in place. The variable parameters 
(i.e. couplings ^ic and poles so,i) are then adjusted to fit the VPI-KSU partial wave T-matrix elements. 



C. Resonance parameter extraction 

Following reference Igl a resonance position is identified with a pole of the scattering T matrix in the complex 
energy-plane. This can only be done for models that can be evaluated for complex values of s, i.e. models that have 
the correct analyticity structure. In the CMB model, the determinant of the H matrix defined by Eq.(15) equals zero 
at the pole position, s = Spoie-, in the complex s-plane. Only the poles located close to the real axis are interpreted 
as resonances. This procedure involves an analytic continuation of (pds) into the complex s-plane for ^m{s) < 0. 
Clearly, the analyticity, defined by the dispersion relation, Eq.(ll), plays a "dynamical" role in finding the resonance 
parameters. This is one of the main differences between our approach and the KSU approach. 

To proceed, we need to evaluate Eqs.(ll) and (12) for complex s. Each T-matrix element has a branch cut beginning 
at the elastic and each inelastic threshold. The branch cut can have a strong effect on the amplitude, in extreme 
cases producing a cusp, e.g. in the TriV — > irN Sn partial wave. Above each threshold, the amplitude is multi-valued. 
This is traditionally described by a Riemann sheet structure pi. The amplitude is continuous when analytically 
continuing to the appropriate new sheet as the value of s crosses the branch cut, but discontinuous when staying on 
the same sheet. The function, 0c(s) defined by Eqs.(9) and (10) is the channel propagator for the "first sheet" of the 
complex s-plane, labeled (j)c,i- At s > Sth.c, it has a discontinuity in its imaginary part determined by unitarity as s 
crosses the real axis. 

The resonance pole is on the "second sheet" in which (j)c{s) = (t>c,ii{s)', it has the same discontinuity as the first 
sheet except for the opposite sign, 



PcJI 



(s + ie) - (jjcjii-s -ie) = (j)cj{s -ie) - (j)cj{s + ie) ior s > Sth (16) 



It is (j)cji{s) which is used in the search of the resonance pole positions. Using the same strategy as CMB p, we 
search the TriV -^ ttN T matrix for poles on the sheet most directly reached from the physical region. As discussed by 
Cutkosky and Wang p3] and elsewhere in the literature, each resonance has additional poles on other Riemann sheets 
associated with each inelastic threshold. However, the pole closest to the physical region is most closely associated 
with the physical characteristics of the resonance. The formalism presented here can be extended to search for poles 
on other sheets and try to distinguish between resonance poles and poles that arise from bound states of composite 
particles. 

A resonance pole is found by searching for a zero in the determinant of the H matrix defined in Eq. (15). Once a 
pole is found, H{s) is diagonalized at the pole position Spoie- This is done to eliminate resonance-resonance interference 
effects when multiple resonances are present in a partial wave. By using the resulting eigenfunctions Xi, the T-matrix 
in the vicinity of the pole can be written as H 



^'^'^^2^+1^ D{7) (^^^ 

ef ^ ' 

The 77c describe the coupling of the resonance to channel c, as defined in Eq. O. This form can be shown to be 
equivalent to the full CMB model. Eq. (|l^ is a general form for a Breit-Wigner resonance shifted by nonresonant 
(background) reaction mechanisms. To determine the resonance properties, we look at the T matrix in the vicinity 
of the pole with a simplifying assumption. The form of the background part, Bab, is assumed to be smooth in the 
immediate vicinity of a pole. This allows us to use the denominator, D{s), to define the Breit-Wigner resonance 
parameters at the pole, ignoring the background. The denominator of the above expression is then matched to a 
relativistic Breit-Wigner form. The full denominator can be written as 

D{s) =r-s- i!Y,yc(pc (18) 

c 

where 



Vc = hd^ 



E 



j=i 



licXi 



(19) 



The real constants, r and v, are defined by equating D{spoie) = 0. The resonance mass comes from the real part of 
the denominator and the width comes from the imaginary part. 

The general form of the Breit-Wigner denominator for a multi-channel situation is 

DbW = {Wl^s -s)-i Wres ^ Tres,c (20) 



where Wres is the mass of the resonance and Tres,c is the decay width of the resonance to asymptotic state c. Thus, 
shifts in the resonance mass can be identified with the real parts of the denominator and the imaginary part can 
be identified with a sum over the partial decay widths of the resonance to various channels, c. To obtain numerical 
values, we use a linear approximation for the real part of the summation in Eq. nS. 

Using the above definitions, the following qualitative statements about the relationship between the resonance 
parameters and the model parameters can be made. The resonance width is approximately equal to twice the 
imaginary part of the denominator at resonance. Because the self energy term in the denominator has strong energy 
dependence, the physical pole position is shifted from the bare pole by an amount that depends of the value of the 
coupling parameter (7) and both the value and shape of ^e(j){s) in the vicinity of the pole. The mass of the resonance 
is further shifted from the physical pole because the mass is determined on the real axis of the complex energy plane. 
Shifts from the real part of the bare pole to the mass can be positive or negative; the size of the shift depends on 
many factors and can be quite large (more information can be found in section M) . 

Eq. (|l^) has an energy dependent pole shift due to 5Re(/)(s). We wish to remove this energy dependence to make 
a connection between the full T-matrix in Eq. (n^ and a standard relativistic Breit-Wigner shape. Making the 

assumption that the "real" part of the term yjj/c^c is linear close to the resonance pole: 

c 

v^e^yccpc f^ a + l3s (21) 

C 

the T-matrix can be re-written in the form of a relativistic Breit-Wigner resonance, without an energy dependent 
pole shift: 

„ Bab -Sab , V^ Vl+0 yi+P .„„^ 

Jab = :^. \' 2^ / _ \ v^ ^ > 



2i 



In Eq. (p2|), the quantitative expressions for the resonance parameters can then be identified in terms only of 
quantities evaluated at the resonance mass : 

^eD{M^^^) = Defines Resonance Mass -^ Mres 

V — jvj-'''^ej)/(7/i ) Defines Resonance Width (22,) 

Y^ _ y^^mip^ — P pggjigg Branching Fraction into channel c 

2^^ ya^mct>a 

In the above equations, D'{s) is the derivative of D with respect to s, 5ftel?'(s) — —(1 + /3). This term accounts for 
the fact that there is an energy dependence to 3fie(/)(s), which shifts the pole position. 

This formulation is identical to the generalized Breit-Wigner form that is the basis of most fitting and theoretical 
models. However, we only use this form in the immediate vicinity of the pole to determine the resonance parameters. 
The data are still fit with the full model and the pole position is then determined from that fit. This type of prescription 
for obtaining the mass, width and partial widths of a resonance is model dependent. It is important to realize that 
all definitions of baryon resonance parameters are model dependent. For a highly elastic and isolated resonance, 
the model is not very important. However, very few resonances fit this description and we are trying to present a 
formulation that minimizes the model dependence. 

D. Relationship to other models 

The CMB model has a number of features that are not included in commonly used models. Here, we present 
a discussion of formulations similar to standard models by making approximations to the full CMB model. These 



simpler formulas will be used in section tV D to show the corresponding model dependence in the extracted resonance 
properties. 

The full CMB model contains a dispersion relation which guarantees the analyticity. The imaginary part of the 
channel propagator function (Eq. n^) is the rclativisticphase space function and the real part is then calculated from 
a dispersion integral. Models which are not analytic jj] include only the phase space, so we set 9fie(/)(s) = for all s 
to simulate them. 



The full CMB model uses a Dyson equation to allow for conversion into open intermediate states (resonant or 
nonresonant) and open asymptotic channels. The bare propagator (Eq. |g) is "dressed" by all the open intermediate 

j] uses the bare propagator of the CMB model in place 



and asymptotic states. The X-matrix formulation P,E^.E1 
of the dressed propagator as a X-matrix rather than a T-matrix and identifies the resonance properties with its 
parameters. A nonanalytic unitary multichannel i^-matrix using a relativistic Breit-Wigner form [E3[ can then be 
constructed for the contribution of resonance R to the reaction between initial state i and final state j: 

\/^rnxJH{s) ^ji^ jRj a/9to(/)j (s) 

sr 

where '^rn(j)i{s) is the product of the form factor and phase space for channel i as defined in Eq. (|l l|) . Mr and Tr^i 
are the mass and partial width for decay to channel i of resonance R; the total width is Tr. These three physical 
quantities are defined by 



j^R(g) ^ v^"'^y^^,^^; iH,i IH,] y ^"^vn^J f24) 



Mr^^^ 



Tfl, 



- M„ (25) 



For the nonrelativistic Breit-Wigner case, 



where Tr is defined by 



r^R( . _ v/Sto(/),(s) 7fl,, 7flj ^5to0j (s) 



rR = 2Y,^rnMsR){lR.Af (27) 



In either case, the corresponding T-matrix can then be found through the standard definition: 

T^K^{I~iK)-^ (28) 

Since both K and T are matrices, there is no simple closed expression for T corresponding to Eqs. (24) and (26) 
except in the single channel case. 

Nonresonant ii'-matrices must also be defined. To maintain unitarity, the full iiT-inatrix is obtained by adding 
all the resonant and nonresonant if -matrices (e.g. [p8[). The corresponding T-matrix includes effects of resonance 
interference and coupling to nonresonant processes, but only on-shell. 

Various theoretical schemes have been built on the K matrix method. At low energies, the characteristics of the 
^ (P33(1232)) and the Sii(1535) have been determined with an effective Lagrangian method |Q. More complete 
formulations have been developed p^,^]. In these models, the mesons and baryons are each fundamental particles. 
Although the number of parameters is reduced from what is required for the more empirical models, there are still 
a number of ambiguities in the construction of the Lagrangian and in the proper development of a multi-resonance, 
multi-channel model. 

III. DATABASE 

We devote a separate section to the database because it plays a critical role in the results we obtain. There are 
many possible asymptotic states that can couple to each resonance. There are also a number of resonances which 
couple weakly to the ttN channel so that the state is only seen in the inelastic data. Just as it is important to include 
various inelastic channels with proper threshold effects in the theoretical model, it is also important to include as 
much of the relevant data as possible with appropriate error bars. 

Although it is best to use the original data, various partial wave analyses producing T matrix representations of the 
data are available. This kind of analysis can be done with much less model dependence than is found in the analysis 
used for determination of baryon resonance properties. Nevertheless, these analyses make choices in the data used in 
the fits and their absolute normalizations since not all data sets are consistent with each other; these choices can add 
error beyond what was in the original data, thus adding uncertainty to the fit. On the other hand, fitting partial wave 
amplitudes allows a simpler fitting strategy- separate fits can be made for each partial wave and less computer time is 
required. This is the same procedure chosen by KSU [Q. We choose to fit the single-energy partial wave amplitudes 
of the VPI group for elastic scattering ||] and the isobar model fits of KSU-VPI for inelastic pion production [g5| . 
We also make a separate partial wave analysis of the ttN -^ rjN data. 



A. ttN elastic data 

Data in this channel is easiest to measure. Therefore, the data are more complete and of higher quality than in the 
other channels. In many resonance parameter analyses, these data have a dominant role in the results. This takes 
advantage of using the best data. However, the inelastic data must be included to find resonances that would not be 
seen in ttN elastic scattering. The Constituent Quark Model l^,^ predicts the existence of a large number of states 
(roughly a number equal to the number of states seen to date) and finds that many of them couple weakly to the ttN 
channel. 

There have been numerous analyses of the elastic data. Since complete experimental results are not yet available, 
theoretical constraints must be employed to get unambiguous fits. Older analyses of CMB ||2^| and KH had 
strong reliance on dispersion relations to generate unique fits. KH80 used fixed-t and fixed-0 dispersion relations 
at many angles. The real parts of these results were later compared with partial- wave fixed-t dispersion relation 
predictions by Koch. ||2l| Although KH80 results are quite noisy close to threshold in the high L partial waves, there 
is good qualitative agreement with the additional constraints. CMB80 uses hyperbolic constraints in the Mandclstam 
variables. 

Although the older works have the best theoretical underpinnings, we use the latest ttN elastic partial wave 
amplitudes of VPI M . The VPI analysis uses a significantly larger data set than was available for the earlier CMB80 
and KH80 elastic analyses. Most of the new points are at W <1600 MeV, but there are also many new results at 
higher energies H]. Thus, a number of older data points with significantly larger estimated errors could be dropped 
from the fit, resulting in an improved fit. These more recent results are consistent with fixed-t dispersion relations 
at 6* = and low W. Although improvements in the VPI analysis are expected, we use this somewhat debatable p^] 
approach in this work since this provides a way to include all of the data published since the older work. Desire for 
an update of the 1980's work has been expressed at conferences for several years, but no such work is in progress. 

B. TvN ^ rnvN data 

For the best available representation of inelastic data, the choices are much more limited than for elastic data. 
Since there are no model-independent methods for the 3-body final states, isobar models are employed. Since no 
recent interpretation of the nN —> t:t:N data is presently available, we use the quasi-two-body channel decomposition 
of Manley, Arndt, Goradia, and Teplitz |^. That work fit an isobar model to the ttN -^ nnN data, isolating the 
contributions of ttA, pN, aN (with a representing the tttt strength in an isoscalar s wave), and 7rA^*(I440) channels. 
Although the aN channel will absorb some of the nonresonant tzttN strength, these choices ignore some of the non- 
resonant ttttN strength and states such as 7rA^*(I520) that might be expected to share strength. They used a data 
sample of about 241,000 events spread over 1320-1910 MeV in W; statistical accuracy is much poorer than for the 
elastic channel and there is no data at the highest energies where resonances are found. All data published after this 
analysis were very close to threshold. It will be clear from our analysis that these data need augmentation and that 
the isobar analysis should be repeated. 

As with the choice in elastic data amplitude, this choice is a compromise between including all the existing data 
in a simple form and using a more proper analysis. An analysis the TriV inelastic data without the assumption of the 
isobar model would be extremely difficult and unjustified with the quality of the present data set. 

In the CMB analysis Iq,^, the ttN -^ nnN data were weighted by a factor of 1/3 smaller than the results of 
Manley, et al. [p5| because they felt the errors were understated. We have a similar attitude and weight the inelastic 
data by a factor of 1/2. Otherwise, the fits to the higher quality elastic data are degraded. 

C. Description of the tvN -^ rjN Analysis 

The most significant inelastic channel at low W is the rjN state; data involving this channel is thought to be crucial 
in analyzing Sii(1535) because the s-wave nN — > rjN cross section is large and rapidly changing close to the resonance 
mass. Unfortunately, the data in this channel is both limited and of uncertain quality. |2g|] The data set with the 
most points close to threshold was published by Brown, et al. [E8|. Clajus and Nefkens p9[ argue that these data 
have unknown errors in the assigned values of W, making them unusable. As mentioned earlier, Batinic et al. |27| ] 
put a very small weight on these data points in their partial wave fit to the nN — > r]N data. To produce partial wave 
amplitudes for this reaction, we reproduce the Batinic et al. fit to the nN -^ ijN data, effectively leaving out the 
controversial Brown et al. ||2ll data set. 



Since there is very little data available for the reaction ttN — > 77A^, we use a simplified version of the full CMB 
model to simultaneously fit the TriV elastic T matrices and the ttN — > rjN data to provide a parameterization for each 
of the partial waves with L <4 contributing to the 77 production cross section. The partial wave analysis followed 
the procedures used by reference |2^] , but used the newest VPI TriV elastic partial wave amplitudes [|| instead of the 
older Karlsriihe Helsinki (KH80) amplitudes. The channels used in the analysis are ttN , rjN, and a dummy channel 
consisting of a fictitious meson with mass chosen so that the dummy channel opens at about the energy where ttN 
inelasticity due to channels such as ttttN start. The mass of that fictitious meson changes from partial wave to partial 
wave. The values used in the analysis are given in Ref. (2^. The partial waves used are the 1=1/2 partial waves 
through G17. All partial wave parameters were varied simultaneously. More details on this procedure given in this 
section can be found in the references [2^,^ . 

For this process, the Sii(1535) resonance makes up most of the total cross section. Therefore, the Sn partial wave 
is the most accurately extracted. The other partial wave amplitudes are smaller and less accurately determined. The 
results of this fit are T-matrices that best model the partial wave data. The best fit is shown as an error band in 
the figures of the results section. For the final fits, 40 data points between threshold and 2.3 GeV were used for each 
partial wave. This insures that these data provide the appropriate contribution to the total chisquare. 

IV. ILLUSTRATIVE EXAMPLES- P33 AND Sn PARTIAL WAVES 

The purpose of this section is to introduce features of the model through the examination of specific partial waves. 
The P33 and Sn partial waves are chosen for this purpose. The P33(1232) or A is the best example of an isolated 
and elastic, i.e. simple, resonance; its characteristics are largely well-established. On the other hand, the Sn partial 
wave has some of the most interesting structure of all the partial waves contributing to TriV scattering. Because of 
this structure, extra care must be taken when extracting resonance parameters for this partial wave. Some of the 
interesting features exhibited in this partial wave are listed below. 

1. There are 2 PDG 4* resonances (Sii(1535) and Sii(1650)) which overlap significantly. 

2. The Sii(1535) has strong coupling to both the TriV and rjN channels and is very near the rjN threshold («1487 
MeV). This produces a strong cusp in the TriV elastic Sn T-matrix element. 

3. There are 7 decay channels that have measurable coupling to the Sn states, each of which has different phase 
space which can cause structure in all of the other channels via unitarity. 

All of the above features are adequately handled in the CMB Model. 

This section will help in understanding the unitarity, analyticity, and other properties of the model, which are 
especially important in the Sn partial wave. A discussion of the cusp associated with the Sn partial wave will also 
be given. Furthermore, this section will present a systematic study of the model dependence of resonance parameters 
extraction by examining results obtained when leaving out various features of the full CMB model. 

A. Single Channel-Single Resonance Case 

The equations from section Q give expressions for a reaction T-matrix element for any number of open scattering 
channels or asymptotic states (rriV, r/N, etc.) and any number of intermediate states (Sii(1535), P33(1232), etc.). 
The T-matrix elements for one asymptotic state and one intermediate state are much simpler. The resulting equations 
can then be apphed with good success to the P33 partial wave near the A(1232) resonance. Despite the complexity 
of the CMB formulation, the T-matrices for an isolated, single channel resonance have a form similar to that found 
in the Breit-Wigner shapes commonly used. 

The one asymptotic state-one intermediate state case requires two parameters (using A to label the P33(1232) 
intermediate state and TriV as the asymptotic state), one coupling ja,ttN and one bare pole energy, sa- Since there 
is only one intermediate state the matrix equations are reduced to scalar equations and the math is simple. There 
is only one channel propagator cjJj^n which is determined analytically with a contour integral. (There is a functional 
form for the channel propagator only in the S-wave case. The results for the TriV channel in an S-wave are shown in 
figure pi) There is also one self energy term according to Eq. 0. 

Sa,a(s) = lA,7vN(l^7rN{s) (29) 

Since the function varies with energy, E does also. The H matrix is then defined by: 



H = (sA - s) - Sa, 



A 



/ \ 2 1, (30) 

= (SA - S) - 7A,,riV<P^W ^ ^ 

The G matrix (the dressed propagator) is then calculated from H as defined by Eqs. (|lj) and (15) in section |l|: 

G=- -^^ -— (31) 

(SA - s) - Ta^-^nVt^n 

TttN,itN is defined using Eq. (n^) and the above equations: 

7A.^JV^"^'/'^JV ,„„. 

(sA - s) - li^^i^4>T,N 

The results of a model calculation for the T matrix is shown in Fig. H with a solid dot at the physical resonance 
mass. The equations produce the characteristic shape of a resonance with a maximum in the imaginary part and a 
zero in the real part at the resonance mass. A similar signal should be seen for strong resonances, but there will be 
shifts in the real data when underlying backgrounds of varying smoothness are included. 

Eq.p3 has the usual form of a relativistic Breit-Wigner resonance as expected, but there are some important 
new features. The real part of the pole is shifted from the bare pole energy (sa) by 7^3ffe(^(s) and the pole gains 
an imaginary part, 7^3m(/)(s) due to couplings of the resonance to the asymptotic state as it propagates. In this 
model, these shifts are energy dependent and come from both unitarity and analyticity requirements. The analyticity 
condition, Eq. (11) then allows an analytical continuation of the above expressions to the complex s plane where the 
physical pole position, mass, and width of the resonance must then be determined by a search discussed in section ||. 

In other contexts, unitarity requirements are satisfied through inclusion of final state interactions. These methods 
add terms to the denominator similar to what comes from the CMB model. However, these models are not analytic. 

Near the resonance peak, Eq. ^ can be expressed in terms of a mass and width identical to a generalized Breit- 
Wigner shape. Thus, the features of the CMB model can be absorbed into effective constants for the case of an 
isolated resonance. 

B. Two Channel- Two Resonance Case 

The equations get rapidly more complicated as the number of resonances and open channels in a partial wave 
increases. The two channel (or asymptotic states) two resonance (or intermediate states) case is still instructive. The 
relevant equations are given below for the L = isospin 1/2 (Sn) partial wave. In reality, this partial wave has two 
strong channels and two strong resonances. The results shown below use realistic parameters for these states, but are 
not meant to be an accurate representation of data because non-resonant processes and other channels are ignored. 
For clarity, channel labels nN and rjN and resonance labels of 1535 (referring to Sii(1535)) and 1650 (referring 
to Sii(1650)) are used. Parameters that must be determined include four coupling strengths (71535, ttAt, 71535, r^w, 
7i650,7rAr, and 7i650,j)w) and two bare poles (S1535 and Si65o)- 

The terms in the S self energy matrix are constructed from Eq. (7): 

Sl535,1535 = 7f535,7rAr'/'''-W + iI535,vN'^vN 

5^1535,1650 = 71535, 7rJV7l650, 7i-Ar</'7rJV + 71535, 7;JV7l65O,?;JV0?)Ar (33) 

Sl650,1650 = 7m50,-^N^-^N + ll^^o^nN^vN 

The S matrix is symmetric by construction so it takes on the simple form: 

S](S) = ( ^1535.1535 Si535,1650 \ (34) 

\ ^1535,1650 ^1650,1650 / 

Fig. shows the S1535, 1535(5) function. It is a weighted sum of ^jrjv and 0,,Af, shown in Fig. 0, and thus carries 
the threshold behavior of both the r]N and irN channels. The other elements have similar qualitative behavior, but 
varying weights of (j) because the two resonances couple to the two asymptotic channels with different strength. 
The H matrix becomes: 

H{s)^ (s,-s)-S(s) 

S1535 — S — Sl535,1535 —5^1535,1650 ^ (35) 

— 5]i535,i650 S1650 — S — Sl650,1650 



The dressed propagator, the G matrix, has a simple form since it is a 2 x 2 matrix: 

G= H-^ 

1 ( 51650 — S — El650,1650 '^1535AQ5a \ (36) 



1^1 \ 5]i535^i650 S1535 — S — Sl535,1535 

where \H\ is the determinant of H: 

\H\ — (si535 — S— Sl535,1535) (si650 — S— Sl650,165o) 



(37) 
The T matrix elements are then defined in section || in Eq. (|^). Explicitly the ttN elastic T- matrix element is: 



^1535,1650 



-'-nN,iTN — ^y^ZTTTZ, '~ 

S1535— S — 2l535,1535 



1650-»-Sl650,1650 
^ 



(38) 



«1650 — S— Sl650,1650- 



-s-E 



1535,1535 



277riV,153577riV,165oSl535,165o3'"l07riV 



(si535—« — Sl535, 1535 )(si650-'5-Si650, 1650 ) — Sj^535_j^ggQ 

Note that all the functions used in this expression depend on s. The above expression is the sum of two resonances 
with an interference term. (If there is only one resonance, this formula simplifies to the same as Eq. (32).) The terms 
that look like single resonances are more complicated than in the previous example. The self energy terms have the 
necessary analytic cuts from the rjN as well as the irN channels and also make the amplitude unitary. There are 
analogous T-matrix elements for the processes: irN -^ rjN and 7]N -^ rjN. 

In this case, two poles must be found in the complex energy plane. Since the two resonances have significant 
interference, the non-diagonal elements of the G matrix are large. This means the second resonance has a contribution 



to the propagator of the first resonance at the pole of the first resonance. Referring to Eq. (38), there are shifts to 
the mass and width of each resonance due to the presence of the other in addition to those due to the asymptotic 
channel couplings which were not seen in the one channel-one resonance case. Here, the H matrix (the inverse of G or 
the "denominator" matrix) must then be diagonalized at each pole to isolate the contributions from each individual 
resonance. 

The ttN elastic and TriV — > rjN T matrices in this partial wave corresponding to the above equation for representative 
parameters (see Table H and O for a complete listing of the relevant values) are shown in Fig. pi The eta threshold 
has a significant effect on the observables. Therefore, the ttN elastic T-matrix element peaks at the r]N threshold 
rather than at the peak of the Sii(1535) resonance. The physical masses are shown as solid dots which are near the 
peak of the cross section and the imaginary part of T for the higher state, but above those positions for the lower 
state. For all channels, the peak of the T-matrix is shifted by a few dozen MeV from the bare mass by the self energy 
terms. (The Sii(1535) resonance shape is also modified by the presence of a strong inelastic threshold.) There is 
also significant interference between the two resonances. Both cross section bumps are unusually narrow compared 
to the physical widths because of interference effects. We find that the Sii(1535) properties are significantly altered 
by the interference with Sii(1650), similar to the findings of Sauermann, Friman, and Norenberg ||33[| in a iiT-matrix 
photoproduction calculation. 

C. Cusp structure in the Sn partial wave 

There are two main causes of the strong cusp structure observed in -kN elastic scattering differential cross section, 
as well as the Sn partial wave amplitudes. The first is that the Sii(1535) resonance couples strongly to both the ttTV 
and rjN channels and the rjN channel threshold (« 1487 MeV) is just below the 1535 pole. The second is that the 
orbital angular momentum is zero [tiN S-wave); thus, the cross section in this partial wave increases linearly with 
momentum. Therefore, through analyticity and unitarity, a cusp structure appears in the -kN elastic channel. 

Figure ra shows the self energy term labeled 21535,1535 for the 2 channel 2 resonance model of Sec. IV B. This self 



energy has the appropriate analytic phase space factors and therefore shows a cusp. Since E enters directly into the 
T-matrix elements, the cusp shows up there as well. Figure H shows the T-matrix elements for elastic scattering 
and p production. Although the cusp structure is apparent in the T matrices for all channels other than that of eta 
production, it is most evident in the ttN elastic channel because most of the decay width of the Sii(1535) is split 
roughly equally between the TriV and r]N channels. All other inelastic channel openings have the potential of creating 
a cusp, but this is by far the strongest case. At present, the ttN data is not of high enough quality to see any other 
cusp structures. 



D. Model Dependence in Analysis of the Sn partial wave 

The variation between this and other models can best be seen in an analysis of the data in the TriV Sn channel, as 
discussed in previous sections. Because of the peculiarities of this channel, data for both TriV — + ttN and TriV -^ rjN 
reactions would be required for a high quality determination of the Sii(1535) parameters and data for TriV -^ ttttN 
would be required for a good determination of properties for the higher energy states. Presently, the inelastic data is 
of much less quality than for the elastic channel. 

To test model dependence in the resonance parameters, we have analyzed all available and subsets of the data for 
the Sii channel with the full CMB model and various approximations that simulate the models employed by various 
other groups. Although the final fits were of similar quality for all cases, there can be large differences in the extracted 
resonance parameters. 

In Tables l^andM, masses, widths and branching fractions are given for the two of the three Sn resonances used 
in this analysis for a number of different model types using various data sets to constrain the fits. The four columns 
on the left of the table describe features of the model used and what data was used in the fit. The five columns on 
the right give the results of the fit for the resonance parameters- mass, width, and branching fractions. 

The Unitarity column labels how unitarity was imposed in the fits. K-matrix means that unitarity was imposed 
by commonly used i^-matrix methods of Moorhouse, Rosenfeld, and Oberlack [E3. For the results shown in the 
tables, we recreate the iiT-matrix fit of Ref. |23] with the appropriate phase space factors (<^(s)) from our model as 



discussed in Sec. [ID. Dyson equation means that unitarity was imposed using the Dyson-equation approach of the 
CMB model. Channels can interact any number of times in various forms before finally decaying because the Dyson 
equation is iterative. 

The column labeled Disp Rel refers to whether or not the dispersion relation was used to make the phase space 
factors in the self-energy terms analytic functions or not. If the dispersion relation is evaluated, the phase space factors 
and hence the self-energies are analytic functions of the square of the CM energy, s; otherwise they are not. When 
using the Dyson equation and not evaluating the dispersion relation (in the fits, we set SRe^ = 0.0), all intermediate 
interactions in the scattering process occur on-shell. Note that this does not affect unitarity. If unitarity is achieved 
through iiT-matrix methods, the mass and widths are direct parameters in the fit. 

The column Res Type describes the form of resonance used. NRBW leieis to a non-relativistic Breit-Wigner shape 

2 

and RBW means a relativistic Breit-Wigner shape was used. For the NRBW case, a resonance has a — ^-r form 

2 

and the 7 couplings have units of ^/energy. For the RBW case, a resonance has a _^_^ -^ form; the 7 couplings 
have units of energy. 

Finally, the column labeled Channels in fit describes the types of data used in the fit. For the ttN case, only the 
VPI TriV elastic Sn T-matrix elements were used in the fit. For the TrA*", 77 A'' cases both the VPI elastic data as 
well as constraints from a partial wave analysis of nN -^ rjN done by this group were used. For cases labeled All, 
ttN -^ {ttN, r/N, ttttN} are all included in the fit. The irirN channel is composed of quasi two body channels {piN)s, 
{P3N)d, (7rA)£), aN, and 7rAf*(1440). Partial wave quasi-two-body T-matrix elements of Manley et al. p3] are used 
for these channels. 

The last row in each table contains results for our full model. All channels are used in the fit, unitarity comes from 
use of the Dyson equation, analytic phase space factors are used, and the bare resonance has a _Ri?W dependence. 

Significant model dependence is seen along with sensitivity to the data used. The lower state (Sii(1535)) has 
significant model dependence as noted above while Sii(1650) tends to have poorer quality data, i.e. missing data or 
data with error bars that are too small. We do not show any results for the third resonance (Sii(2090)) because the 
data quality dominates the fitting, causing wide variation in fit results, e.g. masses vary between 1509 and 2028 MeV. 



The data quality for this partial wave is discussed in detail in section V C . 

We note that without the dispersion relation, the T matrices for the X-niatrix model and the model using the Dyson 
equation for the resonance propagator are equivalent. Therefore, only the AT-matrix results without the dispersion 
relation are given in the table. In general, the choice of relativistic vs. nonrelativistic shape for the bare Breit-Wigner 
resonance does not have a strong influence for an isolated resonance such as the 1535 MeV state. However, the 
1650 MeV state has a weaker signal (in part because of poorer data quality) and the two shapes can produce larger 
differences. Even there, the agreement in the case where all data is used (line 3 for NRBW vs. line 6 for RBW) is 
very good. 

More important differences are found when comparing A'-matrix vs. Dyson equation results with the dispersion 
relation included (e.g. line 6 vs. line 9). The former is close to the model employed by Manley and Saleski M. These 
two models have differences of about 10% in the total width and up to 50% in the branching fractions. 

The most important deviation from the full result comes from the use of a truncated data set. For the 1535 MeV 
state, ignoring the interference with the r]N final state causes the model to fit the Breit-Wigner shape to the cusp 
at the r]N threshold. The VPI work |0] has a very small width for the 1535 MeV state; although the r]N channel 



is mocked up, none of the actual data is used. For even the fuU model, leaving out the ttttN final state data (such 
as was done by Batinic et al. |2^) produces 20% deviations in the branching fractions. We reproduce the updated 
results of the Batinic et al. paper p^ . 

Both the mass and the width of the Sii(1535) tend to increase as more data channels are added into the fit. It is 
interesting to note that the small widths (i.e. < 100 MeV) are in situations where only the ttN elastic data is used 
in the fit. Also, the branching fraction of rjN tends to be smaller than that of TriV unless the analyticity is taken into 
account in the phase space factor. In other words, when the cusp is handled appropriately, the rjN channel becomes 
the dominant decay mode. The Manley-Saleski [Q analysis, which does not address the analyticity issue, is unable to 
match the cusp well and concludes that -kN is the dominant decay mode. 

The Sii(1650) and the third Sn resonance have more random shifts in their resonance parameters. At energies in 
the region of the excitation of the third Sn state, only elastic and 77 production data exist and the signal is not strong. 
Therefore, the fit parameters for the third Sn state are largely determined by fitting non-statistical fluctuations of 
the data. The fitted masses vary widely from case to case and we feel the results with the different models do not 
give information about features of the models. 

E. Elastic Data Dependence 

A major question is the balance between model dependence and data dependence. The previous section showed 
the model dependence using various subsets of the data used in this analysis. In this section, we present results 
for fits of the Sn channel using two different sets of elastic partial wave amplitudes. In Table W\ we compare our 
standard results with the results of a fit using the combination of CMB80 and KII80 data sets. This fit uses the 
same elastic data as Manley and Saleski m except for the additional ttN -^ rjN data. The differences seen in the 
table are somewhat similar (in total width) to but somewhat smaller (in branching fraction) than those presented in 
the previous section, verifying that the major difference between the present work and Manley and Saleski is likely 
due to model dependence. The changes in the results are larger for the second state than for the lowest state, but 
the differences are significant for both states. The differences are much larger for the third Sn, but are not shown 
because of the weak evidence for this state. 

The balance between the branching fraction into TriV and rjN for the lower state are very similar in the two cases, 
but the total width found using the older data is 22% larger. This gives more evidence for the findings of the previous 
section that differences in branching fractions are a primary result of the difference between this model and that used 
by Manley and Saleski. 

F. The nonresonant amplitude 

A major problem in any extraction of resonance parameters is in the careful separation of resonant and nonresonant 
mechanisms. Ideally, an independent calculation would be used for the nonresonant diagrams. Here, we choose a more 
objective strategy and use a smooth background that can then be "dressed" to contain the correct threshold behavior. 
To do this, we add two subthreshold resonances (one repulsive and one attractive) and one very high energy resonance 
in each partial wave as 'bare' propagators. These propagators are then dressed identically to the true resonances. 



More details are given in section II B. 

The separation into resonant and nonresonant components is shown in Fig. 0for various reactions in the Sn partial 
wave. We show the magnitude of the T-matrix for four different final states calculated with final fit parameters. The 
three lines shown correspond to including only nonresonant couplings, only resonant couplings, and all couplings. 
Since the nonresonant and resonant processes are intermixed in the Dyson equation, there is no way to sum them to 
get the full result. 

While some reactions are dominated by resonant processes (e.g. elastic scattering), others are dominated by the 
nonresonant processes (e.g. TriV ~* piN). The resonance excitation must be sampled through a variety of channels 
to provide the full picture. At very high energies (W~1.9 GeV), the lack of data allows the nonresonant processes to 
dominate. Resonance extraction at these high energies is very unclear with the data presently available. 

V. RESULTS AND DISCUSSION 

We have applied the CMB model to the database presented in section p^- TriV single-energy elastic T matrices of 
VPI P], the inelastic T matrices of Manley, Arndt, Goradia, and Teplitz |p5|, and our own partial wave analysis of 



the vrA^ -^ rjN data. The ttN -^ ttttN raw data was not available in time for the present analysis. A reanalysis of 
the ttN — > ttttN data is in progress [ p4| by our group; a more complete analysis can be presented when that work is 
finished. 

The analysis presented here contains features of CMB Q and KSU Q| since we use the formalism of the former and 
a data set similar to that used by the latter. However, the present analysis goes beyond any previously published. We 
present general results and a detailed discussion of the D15, D13, and Sn partial waves. D15 is an excellent example 
of an isolated resonance, but has strong inelastic couplings. The D13 and Sn partial waves each have a strong state 
(well understood for the former and poorly understood for the latter) along with less well understood states; the 
interpretation of most of these states are sensitive to the features of the present model. 

We will compare our results to previous work with an emphasis on works that treat many channels. We also 
compare to the composite results of PDG Q| . 

A. Details of Fitting 

Although the ttN -^ ttttN data set used here is identical to that used by KSU, they do not include the ttN -^ r]N 
data and use older elastic analyses. KSU includes an rjN channel in the Sn partial wave fit, using the requirements 
of unitarity to fix its coupling strength to the resonances. For the elastic channel, they use CMB and KH80 elastic 
amplitudes simultaneously. 

For the inelastic channels involving two pions, we define channels identical to those chosen in Ref. B. They are 



listed in Table III. Each channel has a specific orbital and spin angular momentum; the nomenclature used is given 
in the table. These nN -^ ttttN channels are almost identical to those chosen in Ref. |2^. The only difference is 
in the choice of the mass and width of the fictitious isoscalar-spin meson, taken as 1 GeV for each parameter in 
Ref. |g5[; like Ref. Q, we choose a mass and width of 800 MeV. We agree with their conclusion that the results are 
not sensitive to this choice. 

The number of states sought in each partial wave was the same as used by KSU Q . We also included the same 
number of open channels in each partial as KSU because of the choices in the amplitudes fitted in Ref. (25| . If no 



amplitudes are provided by Ref. |25| for a partial wave, e.g. for F17, an appropriate dummy channel is used to absorb 
the flux. Table VIII labels rows of this type as 'Flux'. 

The number of parameters depends on the number of resonances to be fit and the number of open channels included 
in the fit. For example, the Sn partial wave has three resonances and we fit values for each bare pole and its couplings 
to as many as 8 channels. There are also two subthreshold and one high energy "states" used to simulate background. 
The masses of these states are constrained to be far away from any of the actual resonances. Three parameters, a pole 
and coupling strengths to ttN and rjN are fit for each subthreshold background pole. The high energy state is allowed 
to couple to all open channels. Thus, there are 38 parameters fit in this partial wave. For the D15 partial wave, only 
one resonance (with coupling to four channels) and three background poles were fit, a total of 15 real parameters. 

In the elastic channel, the data quality is reasonable at values of W from threshold to about 2.0 GeV. At values of 
W larger than roughly 1.8 GeV, no inelastic data is available for the ttttN final state. The inelastic data has significant 
fiuctuations. If the data are to be represented by a smooth function (assumed in all analyses), the error bars are 
underrepresented. In fact, the Manley et al. paper |E5[ states that only diagonal errors were included in their output. 
Correlations should be significant in the analysis and can only add to the uncertainties. Although the elastic data 
were able to be fit well, the inelastic data were not. Since the elastic data is of much higher quality than the inelastic 
data, the inelastic error bars were weighted by a factor of 2 lower than the elastic error bars in order to ensure a 
reasonable fit to the elastic data. (The original Cutkosky, et al. paper |^,^ used a factor of 3 to weight the inelastic 
data.) For elastic data, values of y^ / datapoint were 1.7 and 1.6 for the Sn and D15 waves, respectively. For the 
inelastic amplitudes, y^ / datapoint values were 9.2 and 22.2. The x^ values are given for the partial wave amplitude 
values without the extra weighting factor. (We don't quote x^ P^r degree of freedom because the parameters are 
shared between the elastic and inelastic fits.) Although we get qualitatively better fits to the elastic data than Manley 
and Saleski B in most cases, similar quality fits are found for the inelastic amplitudes. (No values of x^ Eire given 
by KSU.) However, the shapes of our inelastic T matrices are qualitatively different in many cases. Specific partial 
waves will be discussed in section |VC| . 

For a complicated multi-parameter fit, errors are difficult to determine because correlations can be significant. The 
partial wave data we use as input quotes only diagonal errors. We include error estimates in all extracted quantities 
due to propagation of errors quoted in the partial wave data. In addition, we add contributions determined from 
additional fits where the background parameterization is varied. 

To allow a full error analysis, separate fits were made to data over a limited energy range to isolate single resonances. 
A 2 X 2 X matrix was used to model 2 channels at a time and a number of fits were made for each resonance to 



determine errors on each of the extracted quantities. In one fit, the errors on the resonance mass and width and 
the error on the largest branching ratio were determined. With -kN as one channel and the channel with the largest 
remaining branching fraction (often ttA) for the resonance under consideration as the second channel, the resonance 
mass, width, and branching fraction together with a simple parameterization of the nonresonant amplitude were fit to 
the data. Errors for the branching fraction to the channels with less coupling strength were determined in fits where 
the mass and width of the resonance were fixed at the final fit value of this work and only the branching fraction 
and background were fit. For the K matrix in these cases, the channel for which the branching fraction was being 
determined was one ii'-matrix channel and the -kN channel or the channel with the largest branching fraction was 
chosen as the second K-matrix channel. In all cases, two fits were done with different simple assumptions for the 
background dependence, either flat vs. linear or linear vs. quadratic. The first type was most often used. Two 
components to the error for each parameter under study were determined from each pair of fits, a relative error for 
each fit parameter from the fit with the largest x^ and the absolute difference between the two determinations of the fit 
parameter. The first component is representative of the statistical error in the data points; the second error is due to 
systematic effects between the different background parameterizations. When two reasonable fits were obtained, the 
components were added in quadrature. For prominent states, the fits were easily made. However, most states required 
many trials to find fits of the appropriate quality. Background shapes, channel choices, W ranges, and weighting of 
the inelastic channel in the x^ determination were all varied until two good fits were obtained. All error bars quoted 
in the tables discussed below were determined from these fits. 



B. General Results 



A list of resonances found in this analysis is given in Tables VII and VIII and compared to the results of KSU M , 
those of CMB [g and the latest recommended values given by PDG ||l|. For various states, many other results 
exist. We make direct comparisons only with the previous results that use multi-channel models and provide error 



bars for their determinations. We show masses and widths in Table VIl and branching fractions to various inelastic 



channels in Table VIII. These are the results of the analysis discussed in sect. IIC. Poles corresponding to most of 
the states found by KSU are found in this analysis, although the properties can be different; states for which there is 
no evidence in the present analysis are Pi3(1900), Fi5(2000), and D33(1940), all of which are 1* states according to 
PDG [Q . Although we obtain fit results for many weak states, a detailed study of the validity of the fit for weak states 
was not attempted. Discussion in this section will be limited to PDG 3* and 4* resonances. In Table IX, we show 
the results for pole p ositions of this analysis. The pole positions have less model dependence than the parameters in 



Tables VII and VIII 



The analysis presented here is not identical to any previous analysis. The data base used is similar to KSU, but the 
formalism has significant differences. The formalism is identical to CMB, but the data base used is quite different. 
CMB used the best representation of data at the time- their own elastic amphtude analysis |20[ and the inelastic 
quasi-two-body amplitudes of the SLAC-Berkeley |g^ and Imperial College |3g] . Batinic et al. [27|] use a truncated 
version of the formalism used here and still another data set. 

Strong isolated resonances that have a strong elastic coupling are fit well with all models and results for the 
resonance parameters, such as the Di5(1675) and Fi5(1680) masses and widths, tend to have close agreement between 
previous results and the new results. PDG gives a range of 15 McV for both masses and a range of 40 and 20 MeV 
for the Di5 and F15 full widths and our values are within these intervals. For the elastic branching fraction, PDG 
suggests a range of 10%; our results are inside the range for F15 and just outside it for D15. 

The benefit of the multichannel analysis is readily apparent for states with a very small elastic branching fraction. 
For example, Di3(1700) and Pii(1710) (both PDG 3* resonances) are not seen in the VPI elastic analysis Q because 
there are no strong signs of the resonance in the elastic T matrix. However, there is a strong resonance signal in the 
ttN — > ttA T matrix in each case. 

For the cases where KSU differs significantly from the consensus of previous results, the S3i(1620) mass and the 
Sii(1650) elasticity, the new results tend to agree with the older values. To test for the dependence on the elastic 
data used in the fit, we re-did fits with the same elastic data used by KSU and found results qualitatively similar to 
our final analysis. 

An unusual feature of the CMB analysis is the large value for the Sii(1535) width, about 40% larger than any other 
analysis. This analysis obtains a width in closer agreement with KSU and KH than CMB. We have been unable to 
reproduce the large width. We note that CMB based their fit on a subsidiary analysis of Bhandari and Chao |M with 
a compatible model. Although the ttN -^ rjN data has changed little since then, the elastic data we use is a global 
fit to all data rather than the single data set they used. The VPI SM95 single energy solution is systematically larger 
and somewhat fiatter than the data used in that 1977 work. The most confusing aspect is that Bhandari and Chao 



quote a full width of 139±33 MeV, much more compatible with the present results than with the full CMB results. 

There has been significant interest and controversy in the properties of states in the Pn partial wave. In large part 
this is because the data in this partial wave has always been poor. In the present fit, x^ I datapoint for this partial 
wave is 3.8 for the elastic channel and ranges from 8.5 to 17.6 for the inelastic channels. With the low quality of the 
existing inelastic data, no effort was made to determine the correct number of Pn states. We find three Pn states 
with properties somewhat different than those previously obtained, although all of our values are within the suggested 
ranges of PDG. In our results, Pii(1440) has a mass that just overlaps the PDG window and one of the largest 
widths obtained. In our fit, this state has a significant contribution from nonresonant interactions; that together with 
the low quality of the present data produces a large estimated error for the values given. KSU results for this state 
show a somewhat smaller width and a much smaller error bar. We can only comment that the KSU model handles 
background quite differently than the present model and they used older ttN elastic data. For the Pii(1710), we find 
mass and width values well within the large PDG ranges; however, KSU has a very large width with a large estimated 
error. Since this state sits on the tail of the Pii(1440) and doesn't have a strong signal in any channel, it is clear that 
the properties of the 2 lowest Pn states are closely coupled and multichannel analyses are very appropriate. A later 
analysis of Cutkosky and Wang 1 32 1 of this partial wave using the CMB model compared results obtained with the 
VPI (SM89) and CMB80 partial wave analyses. They find qualitatitively similar results (large width for the Roper 
and small width for the 1700 MeV state) to those obtained here. Evidence for the highest energy Pn state in the 
present data is very poor because only elastic data exists in the appropriate energy range. 

Other unusual cases found in this new analysis include the P3i(1910) and S3i(1900) states. The P3i(1910) is found 
at significantly higher mass and has significantly larger width than previous determinations. Since this state is at 
high mass, the inelastic data is very important in determining its properties. However, there is almost no existing 
inelastic data in this partial wave other than a few points in the ttN -^ 7riV*(1440) reaction. PDG has given this state 
a 4* rating, but with the present data this rating should be downgraded. Although VPI doesn't find the S3i(1900) 
state, it is very prominent in KSU and PDG. We find the mass at 1802 MeV and a width of 48 MeV (with a large 
estimated error) while KSU finds 1950 MeV and 263 MeV. This significant difference is largely due to the elastic data 
sets used. There is a strong bump at about 1900 MeV in the elastic T matrices used by KSU which has vanished in 
the VPI partial wave amplitudes. 

For the branching fractions presented in Table VIII, the only recent result is KSU, which is presumably weighted 
heavily in the PDG listings (also shown in the table). We divide the pN and ttA channels into the appropriate spin 
channels since the component angular momenta can sometimes have more than one value. For D13, the pN spin can 
be 1/2 with orbital angular momentum 2 or 3/2 with orbital angular momentum or 2. Since only the spin 3/2 
orbital angular momentum case was found to be important in the irN -^ ttttN isobar analysis |^ , this is the only 
pN channel we include for this partial wave. The ttA channel can couple with orbital angular momentum of or 2 in 
this partial wave and both possibilities are included in the fits. 

As with KSU, uncertain fits due to underestimated error bars and/or missing data make interpretation difficult in 
some cases. We quote values for this analysis and give estimated errors for each quantity. 

Since we use the same inela stic T matrices as KSU for input to the fit, the results should be qualitatively similar. 
Based on our study in section IV D, we feel there is roughly 20% difference between the KSU and Pitt-ANL results 
due to model dependence in the most sensitive quantities. The poor fits make this true less often than might be 
expected. For Fi5(1680), the elastic branching fraction is ~ 70% in both analyses and the largest inelastic channels 
are ttA P-wave and (7r7r)sA^ in both; agreement is within errors for the largest values. For D33(1700), the elastic 
branching fraction is small and the inelastic strength in concentrated in the ttA s-wave, so we are in close agreement 
with KSU. However, S3i(1620) is a strong state where agreement is not good. There are 2 strong inelastic channels, 
piN and ttA D-wave. Although there is agreement in the strength of piN, the elastic strength is much smaller for 
KSU and the ttA branching fraction is of course larger. The full widths seen are only 15% different and the estimated 
errors overlap. The x^ P^r data point for the S31 partial wave in the elastic channel is 2.1 and ranges from 4.3 to 17.2 
for the inelastic channels. 

The calculation of Capstick and Roberts uses a relativized quark model to calculate resonance masses and a '^Po 
model for creation of qq pairs |18|. They fit the two quark level coupling parameters to the nN decay amplitudes 
of the non-strange resonances rated by PDG as 2* or better and then predict the remaining decay amplitudes for 
various meson+baryon final states (including many final states for which there are no data). Although their primary 
purpose was to look for signatures of "missing" quark model resonances, we note that the qualitative agreement with 
our analysis is satisfactory. There are notable successes such as Di3(1520) where quantitative agreement comes with 
all significant decay channels and notable failures such as Sii(1535) where both TriV and rjN decay widths (and the 
full width) are overestimated by a factor of 4. 



C. Detailed Discussion - D15, D13, and Sn Partial Waves 

Presented in this section are detailed results for three representative partial waves. It includes figures of all the 
channel T matrices and a discussion of the results. These partial waves have reasonabl e qua lity data and contain both 



prominent and less prominent states. Numerical data can be found in Tables |^ 0, |VII] . and \Vl% The T matrices 
are found in figures for the D15 partial wave, M for D13, and O for Sn. In each we show the relevant phase shift 



amplitudes Ig,^ along with the fit of Manley and Saleski |^ (dashed lines) and our final fit (solid line) . Only data 
up to W=2.15 GeV were used in the fits because the data at higher values of W are of diminished quality. Since all 
the data in each partial wave were simultaneously fit with the requirements of unitary, not all channel data shown in 
the figures is fit equally well. In general, the fits to the elastic data are good {x^ / datapoint ~ 2) and the fits to the 
inelastic data are poor {x^ / datapoint > 10) for both the KSU fits and the present results, as discussed above. 

The D15 partial w ave co ntains a strong (PDG H 4*), isolated resonance for which the results of the one resonance 
example (see section |VA) apply fairly well. In this case, both the elastic and ttN -^ irA^Dwave) inelastic channels 



are prominent in the D15 partial wave. Since neither threshold is nearby and there is no interfering resonance, the T 
matrix in each of these channels has the characteristic shape of an isolated resonance with a peak in the clastic and 
inelastic cross section near the resonance mass. The non-resonant amplitude is a smooth function. The KSU and 
present fits are similar in the ability to match features in the data in the two strong channels near the cross section 
peak at W ~ 1680 MeV. The two analyses differ in their description of the high energy data {W > 1.8 GeV), but the 
data are very sparse there. Neither model fits the smaller p production channels well. The p^N channel shows that 
the model expectations (based on the requirements of unitarity and the data in the other channels) don't match the 
data in this channel and a poor fit results. 

Essential ly a ll analyses for the Di5(1680) resonance give similar values for the mass, full width, and elastic fraction 
(see Table VII). Even though VPI pj accounts for the inelasticity by using a dummy channel, the fit parameters 
for this state are similar to those obtained in the present fit. PDG ||l| gives quite small error bars for the mass and 
width of this state, reflecting the unanimity of the fitting results for this state. It is encouraging that the complicated 
features of the present model are not important for the simple case presented in this partial wave. 

The D13 partial wave contains a 4* state at 1520 MeV (PDG standard value, the actual mass is slightly different), 
a 3* state at 1700 MeV, and a 2* state at 2080 MeV. Since the lowest two resonances have significant overlap in 
energy, the features of the CMB model are important for this partial wave. The lowest state is highly elastic, but also 
shows up prominently in {p3N)s and (7rA)5 final states. Note that the peak in the imaginary T matrix is inverted. 
As already mentioned, the second state is barely seen in the elastic channel, thus is not found in the VPI |^ analysis. 
Both KSU and the present analysis find most of the decay strength in the ttttN final states with less than 10% of the 
decay strength to the elastic channel. KSU differs with this model in the distribution of the inelastic strength. While 
KSU finds most of the strength in pN, this model finds it in ttA. The KSU model fits available inelastic data better. 
Feuster and Mosel g] find a systematically lower width than other models for the lowest state, but their width for 
the 2nd state is unusually large. Evidence for the third state is very weak in this analysis; there is only elastic data 
a.t W > 1.9 GeV and there is very little evidence of structure. It is best fit with a very large width. 

The 3 resonances have varying strength in the different channels. As a result, the fitted curve also changes sig- 
nificantly from channel to channel. In the ttA channels, the shapes are complicated because the 1700 MeV state 
is important. The imaginary part of the T matrix has a dip at the 1520 MeV state and a peak at the 1700 state. 
Nevertheless, this model and the KSU model fit the features well and the fitted parameters for the higher energy state 
agree within stated errors. 



The Sii partial wave was already discussed in section IV with regard to the significant model dependence seen. As 
a result of this model dependence, the Sii(1535) full width has been quoted as 66-250 MeV based on different analyses 
of similar data (almost all fits are based prominently on the elastic ttN channel). The relatively small width of the 
Sii(1535) in this analysis (112±30 MeV) is determined in large part by a significant overlap of the Sii(1535) with the 
Sii(1650). This overlap causes a large interference effect. Although a similar effect was seen in one description of the 
7P —> rm data that uses a formalism similar to the Pitt-ANL model |33] , many models based on a T- matrix formalism 
(e.g. I38|4l0|) have a formalism where it is difficult to account for more than one resonance in a partial wave. 

The multichannel iiT-matrix analysis of Feuster and Mosel |5| gets values of the full width in the range of 151-198 
MeV for different form factors. However, their model doesn't handle the rjN cusp as well because it lacks the analytic 
phase space, doesn't model the ttttN inelasticity well because a zero width particle is used, and doesn't include off-shell 
intermediate state scattering effects because of the iiT-matrix approximation. 

Although 8 channels are included in the fit, many of them turn out to have small coupling to Sii(1535) (in agreement 
with PDG j|]). The two major channel couplings for Sii(1535) are nN and rjN. Therefore, another determining factor 
in the full width is the total cross section for the Tr~p — > rjn in Fig. U2 and the corresponding T matrix in Fig. Ob. 
The cross section of Fig. h2 is made up of TriV — > rjN T-matrix amplitudes for Sn through G17, however its main 



component is Sn. Although the shape of the cusp (see Fig. M) does carry information about the rjN channel, the 
overall data quality isn't satisfactory. We await higher quality data for a more conclusive determination of this 
amplitude. 

The branching fraction of the Sn (1535) to the ttN channel is at the low end of the PDG range. A surprising feature 
of the Sii(1535) has been its unusually large decay width to rjN. We find strong evidence for a value at the high 
end of the PDG range. We believe our result because our model fully accounts for the threshold enhancement (cusp 
effect) due to the rjN channel opening and the interference of the 2 overlapping resonances. Since no other model has 
such a complete formulation, they can have model-dependent systematic errors as a result. 

The properties of Sii(1650) are largely determined through the prominent bump in the elastic channel. The elastic 
T matrix has a peak in the imaginary part at about the right energy and the real part is decreasing at that energy; 
a strong resonance signal. The strongest inelastic signal comes in the psN channel. The decay width to tjN is small 
in this analysis, in agreement with previous analyses. The striking difference between branching fractions to rjN has 
been another unusual feature of these states. The existence of the third Sn resonance is very weakly supported by 
the existing data. In our analysis, this state is most often determined by the need to fit non-statistical features in the 
partial wave amplitudes. Further data will be required to sort out the question of whether or not it exists. 

The role of the inelastic rjN channel is very important in the determination of resonance parameters for the Sn (1535) 
state, as shown in some detail in section |^ Unfortunately, the existing data for TriV — > rjN is sparse and of uncertain 
quality p^ (see Fig. |l2|). Recently published data pi for ^p -^ rjN has been interpreted as evidence for a very large 
width (~ 200 MeV). We do not include the Krusche et al. data in this analysis. However, very simple resonance 
models [^,Q were used to determine the large width. Since the large body of t:N data was not fitted with their 
model, possible inconsistencies exist in their parameterization. Although the simple model has validity for the rjN 
coupling to Sii(1535), we feel there will still be significant effects from interference with Sii(1650) and coupled channel 
effects (e.g. -yN -^ S'ii(1650) ^ t:N ^ S'ii(1535) ^ TyiV). 

D. Observables 

T-matrix amplitudes were fit using the Pitt-ANL model. We do not fit directly to data; rather, we fit the partial 
wave amplitudes. It is then interesting to compare the actual experimental data to our calculated observables to 
verify that nothing has been lost by not fitting to the actual data points. 

The / and g spin non-flip and spin flip amplitudes are defined in terms of the Pitt-ANL T-matrix elements for TriV 
elastic scattering by p^ : 

9(0) = lj:ZATi+-Te-)p;icose) ^•'''> 

where Pi and P^ are the Legendre polynomials and their first derivative with respect to cos 9, respectively. The Ti± 
are the amplitudes for a particular ttN scattering charge channel, where the £± refers to the total spin J = £±^, and 
q is the center of mass momentum of the final state pion. The Ti± are defined by Eq. (1^) for the nN — > ttN channel. 
The differential cross section (^) and polarization P are: 
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^ = |/(0)P + |sin%(0)|^ 
§P = 2Qm[fi0)g*{e)]sinl 



(40) 



Figure O shows a plot of the total cross section vs. the Pitt-ANL calculation for the process n^p -^ ■K^n. Figure 14 
shows plots of the differential cross section, ^, at a number of different energies for the process 7r+p -^ n^p. Figure 15 
shows plots of the nucleon recoil polarization, P, for a number of different energies for the process n~p —^ TT~p. In 
each case, the calculation using final fit parameters is in good agreement with the data. For the P data, the error 
bars are sometimes large. 

For the plot in figures |l^, |l^. and |l5| the data is taken from VPI's SAID database. It is a large collection of 
scattering data for a number of different reaction channels, including TriV elastic and charge exchange reactions. 

Based on our studies, we can conclude that our resonance parameters are very close to what would have been 
obtained if the fit had been made to the original data. 

VI. CONCLUSIONS 

We have presented a new study of excitation of baryon resonances through pion- nucleon interactions. An expanded 
version of the Carnegie-Mellon Berkeley (CMB) was used in this work. We added the rjN channel to the calculation 



and used the latest data sets. The model satisfies many of the desirable theoretical constraints- two and three body 
unitarity, time reversal invariance, and analyticity. Because the model is based on a separable interaction, it is linked 
to a Hamiltonian approach, unlike many of the models previously used pUl]. The model has a key feature of allowing 
large numbers of asymptotic channels and more than one resonance per partial wave. Since analyticity is satisfied 
and all inelastic thresholds are included, structure in the physical observables due to nonresonant effects are carefully 
treated. Resonances can then be found as poles in the complex s plane. As in the original CMB model, non-resonant 
processes are included as very wide resonances at very high or very low energies which are then subject to the same 
threshold effects as the resonances. This empirical approach is well established in potential scattering theory, but 
could be replaced later by theoretical calculations based on meson exchange models, e.g. |Q. 

Rather than deal with the complexities of a direct fit to the cross section and polarization data, we use the partial 
wave amplitudes of the VPI and KSU groups |§|,^ and our own fit to the ttN — > r]N data. A re-evaluation of the 
TriV —f ttttN data is in progress by our group. 

The features of this model are expected to be very important in the Sn partial wave where there are two overlapping 
resonances at 1535 and 1650 MeV and a strong threshold opening. The lower state pole and the jyN channel opening 
close together and the quantum mechanical interference of the two states is found to be significant. We show how the 
interpretation of these states has very strong model dependence and how relaxing the features of the present model can 
skew the fitting results for the width and branching fractions to ttN and r/N. We get a new result for the branching 
fraction to rjN at the high end of what has been previously published for Sii(1535), making its characteristics even 
more unusual. 

The general fitting results are presented in Tables VII and VIIl| . We searched for the same states seen in the previous 



analysis of Manley and Saleski ^ . This simplification was made because the data quality wasn't high enough for a 
valid search for weakly excited states. Although values of x^ / datapoint are roughly 2 for elastic fits, they are 10 or 
more for fits to the nN -^ ttttN T matrix data. This is perhaps due to the lack of a treatment of correlations in the 
error analysis of the pion production data. Modern experimental techniques could greatly improve the inelastic data 
set. We strongly encourage new measurements of the inelastic channels. A significant effort was made to determine 
error bars reflecting both the estimated errors in the T matrices used in the fits and the differing choices of background 
energy dependence. 

We find strong and isolated states (e.g. P33(1232) and Di5(1675)) with very similar parameters as previous analy- 
ses. [^,E|,||j^]. States with strong model dependence such as Sii(1535) or with significant changes in the data set such 
as S3i(T900) and Di3(1700) get quite different results. We find a full width of Sii(1535) at the low end of the range 
of previous values, especially with respect to the interpretations of recent jp — > rjp data Ell]. The results presented 
here do not include the photoproduction data; rather, the present results are dependent on a fairly weak data set for 
iT~p —^ rjn. Nevertheless, the large widths obtained in analyses highly dependent on the Mainz eta photoproduction 
data are only loosely coupled to the large data set used in this study. 

We are in the process of further extending the model to include photoproduction and electroproduction reactions. 
The Born terms are included for production of pions and etas and the resonance spectrum is being refit. Results of 
the present paper will then be subject to change. 
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TABLE I. Model parameters for a test calculation in the Sn partial wave. These are the bare mass, the physical mass 
and width, and the dressed pole location for the two resonances included in the calculation. No background contributions are 
included. 



state 


bare 
mass 

(MeV) 


physical 
mass 

(MeV) 


full 
width 

(MeV) 


physical 
pole 

(GeV) 


Sii(1535) 
Sii(1650) 


1500 
1650 


1518 
1680 


83 
200 


1.507-0.0321 
1.676-O.lOli 



TABLE II. Model parameters for test calculation in the Sn partial wave. These are the coupling parameters and branching 
fractions (b.f.) between the two asymptotic channels and the two resonances included in the calculation. No background 
contributions are included. 



state 



(GeV) 



-yvN 
(GeV) 



TriV 
b.f. 
(%) 



b.f. 

(%) 



Sii(1535) 
Sii(1650) 



0.5 
0.9 



0.85 
-0.04 



42 
96 



58 
4 



TABLE III. Channels included in this analysis. Since isospin symmetry is assumed to simplify the analysis, all charge states 
of a particle have the same mass. We choose channels very similar to those used by KSU with identical nomenclature. All 
channels used in fits have fixed spin, orbital, and total angular momentum. In the text, the pN channel is denoted by {psN)t 
with s equal to twice the total spin of the pN system and I giving the orbital angular momentum in spectroscopic notation. 
The orbital angular momentum of the ttA is also ambiguous, so it is also given as a subscript. 



channel 



baryon 
mass 



baryon 
width 



meson 
mass 



meson 
width 



TriV 

■nN 

pN 
<tN 

ujN 
ttA 
ttN* 
KK 



939 
939 
939 
939 
939 
1232 
1440 
1116 



115 
200 



139 
549 
770 
800 
782 
139 
139 
498 



153 
800 



TABLE IV. Model dependence for the Sii(1535). See text for details. 



Unitarity 


Disp. 


Res. 


Channels 


Mass 


Width 


ttN 


'/yN 


ttttN 




Rol. 


Type 


in fit 


MeV 


MeV 


% 


% 


% 


K-Matrix 


NO 


NRBW 


ttN 


1518 


87 


43 


6 


51 


K-Matrix 


NO 


NRBW 


ttN , 77N 


1532 


108 


45 


39 


16 


K-Matrix 


NO 


NRBW 


ALL 


1535 


126 


42 


44 


14 


K-Matrix 


NO 


RBW 


ttN 


1514 


84 


35 





65 


K-Matrix 


NO 


RBW 


ttN , 77N 


1533 


110 


44 


40 


16 


K-Matrix 


NO 


RBW 


ALL 


1534 


125 


42 


43 


15 


Dyson eq. 


YES 


RBW 


ttN 


1531 


72 


16 


62 


22 


Dyson eq. 


YES 


RBW 


ttN , jyN 


1526 


114 


36 


41 


23 


Dyson eq. 


YES 


RBW 


ALL 


1542 


112 


35 


51 


14 



TABLE V. Model dependence for the Sii(1650). See text for details. 



Unitarity 


Disp. 


Res. 


Channels 


Mass 


Width 


ttN 


J7N 


ttttN 




Rol. 


Type 


in fit 


MeV 


MeV 


% 


% 


% 


K-Matrix 


NO 


NRBW 


ttN 


1645 


233 


35 


49 


16 


K-Matrix 


NO 


NRBW 


ttN , ?7N 


1689 


225 


67 


31 


2 


K-Matrix 


NO 


NRBW 


ALL 


1694 


259 


72 


16 


12 


K-Matrix 


NO 


RBW 


ttN 


1682 


161 


78 


5 


17 


K-Matrix 


NO 


RBW 


ttN , ?7N 


1692 


233 


75 


15 


10 


K-Matrix 


NO 


RBW 


ALL 


1690 


229 


65 


25 


10 


Dyson eq. 


YES 


RBW 


ttN 


1692 


138 


65 


35 





Dyson eq. 


YES 


RBW 


ttN , 77N 


1676 


104 


54 


45 


1 


Dyson eq. 


YES 


RBW 


ALL 


1689 


202 


74 


6 


20 



TABLE VI. Sii Resonance parameters from fits to different elastic data sets. All results use the nN — » rjN and nN -^ nnN 
partial wave amplitudes as in the full analysis. Comparisons are made for results using the VPI and a mixture of the CMB80 
and KH80 elastic partial wave amplitudes (as was done in Ref. [4]). The first two lines in the table give results for Sii(1535) 
and the last two lines give results for Sii(1650). Although all channels were used in the fit, only the total TvnN branching 
fraction is given. 



Elastic 


Mass 


Width 


ttN b.f. 


r/N b.f. 


TvnN b.f. 


Data Set 


MeV 


MeV 


% 


% 


% 


VPI 


1542 


112 


35 


51 


14 


CMB/KH80 


1535 


137 


35 


53 


12 


VPI 


1689 


202 


74 


6 


20 


CMB/KH80 


1691 


222 


58 


15 


27 



TABLE VII. Results for masses, full widths, and elastic fractions for all resonances found in this analysis. All resonances 
found in the KSU analysis were searched for, but not all were found. No attempt was made to find new resonances because 
then data quality is not good enough for a new search. See text for more details. 



Resonance 



Mass 
(MeV) 



Width 

(MeV) 



Elasticity 

% 



Reference 



Sii(1535) 



1542(3) 

1534(7) 

1520-1555 

1550(40) 



112(19) 
151(27) 
100-250 
240(80) 



35(8) 
51(5) 
35-55 
50(10) 



Pitt-ANL 
KSU 
PDG 
CMB 



Sii(1650) 



1689(12) 

1659(9) 

1640-1680 

1650(30) 



202(40) 
173(12) 
145-190 
150(40) 



74(2) 
89(7) 
55-90 
65(10) 



Pitt-ANL 
KSU 
PDG 
CMB 



Sii(2090) 
* 



1822(43) 

1928(59) 

^ 2090 

2180(80) 



248(185) 
414(157) 

350(100) 



17(3) 
10(10) 

18(8) 



Pitt-ANL 
KSU 
PDG 
CMB 



Pii(1440) 



1479(80) 
1462(10) 
1430-1470 
1440(30) 



490(120) 
391(34) 
250-450 
340(70) 



72(5) 
69(3) 
60-70 

68(4) 



Pitt-ANL 
KSU 
PDG 
CMB 



Pii(1710) 



1699(65) 

1717(28) 
1680-1740 
1700(50) 



143(100) 

478(226) 

50-250 

90(30) 



27(13) 
9(4) 
10-20 
20(4) 



Pitt-ANL 
KSU 
PDG 
CMB 



Pii(2100) 
* 



2084(93) 

1885(30) 

« 2100 

2125(75) 



1077(643) 
113(44) 

260(100) 



2(5) 
15(6) 

12(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



Pi3(1720) 



1716(112) 
1717(31) 
1650-1750 
1700(50) 



121(39) 
383(179) 
100-200 
125(70) 



5(5) 
13(5) 
10-20 
10(4) 



Pitt-ANL 
KSU 
PDG 
CMB 



Di3(1520) 



1518(3) 

1524(4) 

1515-1530 

1525(10) 



124(4) 
124(8) 
110-135 
120(15) 



63(2) 
59(3) 
50-60 

58(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



Di3(1700) 



1736(33) 
1737(44) 
1650-1750 
1675(25) 



175(133) 

249(218) 

50-150 

90(40) 



4(2) 
1(2) 
5-15 
11(5) 



Pitt-ANL 
KSU 
PDG 
CMB 



Di3(2080) 



2003(18) 

1804(55) 

^ 2080 

2060(80) 



1070(858) 

447(185) 

300(100) 



13(3) 
23(3) 

14(7) 



Pitt-ANL 
KSU 
PDG 
CMB 



Di5(1675) 



1685(4) 

1676(2) 

1670-1685 

1675(10) 



131(10) 
159(7) 
140-180 
160(20) 



35(1) 
47(2) 
40-50 
38(5) 



Pitt-ANL 
KSU 
PDG 
CMB 



Fi5(1680) 



1679(3) 

1684(4) 

1675-1690 

1680(10) 



128(9) 
139(8) 
120-140 
120(10) 



69(2) 
70(3) 
60-70 
62(5) 



Pitt-ANL 
KSU 
PDG 
CMB 



Fi7(1990) 



2311(16) 

2086(28) 

^ 1990 

1970(50) 



205(72) 
535(117) 

350(120) 



22(11) 
6(2) 

6(2) 



Pitt-ANL 
KSU 
PDG 
CMB 



Gi7(2190) 



2168(18) 

2127(9) 

2100-2200 

2200(70) 



453(101) 
547(48) 
350-550 
500(150) 



20(4) 
22(1) 
10-20 
12(6) 



Pitt-ANL 
KSU 
PDG 
CMB 



S3i(1620) 



1617(15) 

1672(7) 

1615-1675 

1620(20) 



143(42) 
154(37) 
120-180 
140(20) 



45(5) 
9(2) 
20-30 
25(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



S3i(1900) 



1802(87) 
1920(24) 
1850-1950 
1890(50) 



48(45) 
263(39) 
140-240 
170(50) 



33(10) 
41(4) 
10-30 
10(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



P3i(1750) 
* 



1721(61) 
1744(36) 
^ 1750 



70(50) 
299(118) 



6(9) 

8(3) 



Pitt-ANL 
KSU 
PDG 



P3i(1910) 



1995(12) 
1882(10) 
1870-1920 
1910(40) 



713(465) 
239(25) 
190-270 
225(50) 



29(21) 
23(8) 
15-30 
19(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



P33(1232) 



1234(5) 

1231(1) 

1230-1234 

1232(3) 



112(18) 
118(4) 
115-125 
120(5) 



100(1) 
100(0) 
98-100 
100(0) 



Pitt-ANL 
KSU 
PDG 
CMB 



P33(1600) 



1687(44) 
1706(10) 
1550-1700 
1600(50) 



493(75) 
430(73) 
250-450 
300(100) 



28(5) 
12(2) 
10-25 

18(4) 



Pitt-ANL 
KSU 
PDG 
CMB 



P33(1920) 



1889(100) 
2014(16) 
1900-1970 
1920(80) 



123(53) 
152(55) 
150-300 
300(100) 



5(4) 
2(2) 
5-20 
20(5) 



Pitt-ANL 
KSU 
PDG 
CMB 



D33(1700) 



1732(23) 
1762(44) 
1670-1770 
1710(30) 



119(70) 
599(248) 
200-400 
280(80) 



5(1) 
14(6) 
10-20 
12(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



D35(1930) 



1932(100) 
1956(22) 
1920-1970 
1940(30) 



316(237) 
526(142) 
250-450 
320(60) 



9(8) 
18(2) 
5-20 
14(4) 



Pitt-ANL 
KSU 
PDG 
CMB 



Das (2350) 
* 



2459(100) 

2171(18) 

^ 2350 

2400(125) 



480(360) 
264(51) 

400(150) 



7(14) 
2(0) 

20(10) 



Pitt-ANL 
KSU 
PDG 
CMB 



F35(1752) 



1724(61) 
1752(32) 



138(68) 
251(93) 



0(1) 
2(1) 



Pitt-ANL 
KSU 



F35(1905) 



1873(77) 
1881(18) 
1870-1920 
1910(30) 



461(111) 
327(51) 
280-440 
400(100) 



9(1) 
12(3) 
5-15 

8(3) 



Pitt-ANL 
KSU 
PDG 
CMB 



F37(1950) 



1936(4.5) 

1945(2) 

1940-1960 

1950(15) 



245(12) 
300(7) 
290-350 
340(50) 



44(1) 
38(1) 
35-40 
39(4) 



Pitt-ANL 
KSU 
PDG 
CMB 



TABLE VIII. Results for decay branch ing ratios of all resonances found in this analysis, 
percentage of the full width found in table VII . 



Fractions are expressed as a 



Resonance 



Channel 



Pitt-ANL 



KSU 



PDG 



Sii(1535) 



ttN 
r?N 
piN 

(/03N)d 

(7rA)o 

(aN)p 

7rN*(1440) 



35(4) 
51(5) 
2(1) 
0(1) 

1(1) 
2(1) 
10(9) 



51(5) 
43(6) 

2(1) 

1(1) 
0(0) 

1(1) 
2(2) 



35-55 

30-55 

0-4 

0-1 
0-3 

0-7 



Sii(1650) 



ttN 
?7N 
PiN 

(P3N)d 

(7rA)o 

(aN)p 

7rN*(1440) 



74(2) 
6(1) 

1(1) 
13(3) 

2(1) 

1(1) 
3(1) 



89(7) 
3(5) 
0(0) 
3(2) 
2(1) 
2(2) 
1(1) 



55-90 
3-10 
4-14 

3-7 
0-4 
0-5 



Sii(2090) 



ttN 
PiN 

(P3N)d 

(7rA)o 

(<tN)p 

7rN*(1440) 



17(7) 
41(4) 
36(1) 
1(1) 
1(1) 
2(1) 
2(1) 



10(10) 
0(3) 
49(22) 
0(1) 
6(14) 
4(10) 
30(22) 



Pii(1440) 



ttN 

?7N 

piN 
(7rA)p 
(<tN)s 

KA 



72(2) 
0(1) 
0(1) 
16(1) 
12(1) 
0(1) 



69(3) 



22(3) 
9(2) 



60-70 
0-8 

20-30 
5-10 



Pii(1710) 



ttN 

?7N 

PiN 
(7rA)p 
(aN)s 

KA 



27(4) 
6(1) 
17(1) 
39(8) 

1(1) 
10(10) 



9(4) 

3(7) 
49(10) 

2(4) 
37(10) 



10-20 

5-25 
15-40 
10-40 
5-25 



Pii(2100) 



ttN 

rjN 

PiN 
(7rA)p 
(aN)s 

KA 



2(1) 
61(61) 

4(1) 

2(1) 

10(1) 

21(20) 



15(6) 

27(79) 

24(18) 

32(71) 

2(6) 



Pi3(1720) 



ttN 

r?N 

piN 

Flux 



5(5) 
4(1) 
91(1) 
0(1) 



13(5) 
87(5) 



10-20 
70-85 



Di3(1520) 



ttN 

(P3N)S 

{nA)s 



63(1) 
0(1) 
9(1) 
11(2) 
15(2) 



59(3) 

21(4) 
15(4) 
5(3) 



50-60 

15-25 
5-12 
10-14 



Di3(1700) 



Di3(2080) 



Di5(1675) 



Fi5(1680) 



Fi7(1990) 



Gi7(2190) 



S3i(1620) 



S3i(1900) 



P3i(1750) 



P3i(1910) 



(<7N)p 

ttN 
ijN 

(P3N)S 

(7rA)c 
{nA)s 
(aN)p 

ttN 

r?N 
(p3N)s 
{nA)D 
(7rA)s 
(aN)p 

ttN 

piN 

(P3N)d 

(7rA)i5 

ttN 
r?N 

(P3N)P 
(p3N)p 

(7rA)p 
(7rA)p 

(<7N)o 

ttN 

r?N 

Flux 

ttN 
r?N 

(P3N)d 

(c<J3N)d 

ttN 
piN 

(P3N)d 

(7rA)o 
7rN*(1440) 

ttN 
piN 

(P3N)d 

(7rA)i, 
7rN*(1440) 

ttN 

7rN*(1440) 

Flux 

ttN 

7rN*(1440) 

Flux 



1(1) 

4(1) 

0(1) 

7(1) 

79(56) 

11(1) 
0(1) 

13(2) 

0(2) 

6(6) 

17(10) 

40(10) 

24(24) 

35(2) 
0(1) 
0(1) 

1(1) 
63(2) 

69(1) 
0(1) 
3(1) 
5(1) 

1(1) 
14(3) 
9(1) 

22(3) 
0(1) 

77(77) 

20(1) 

0(1) 

29(28) 

51(51) 

45(1) 
14(3) 

2(1) 
39(2) 

0(1) 

33(6) 
30(2) 

5(1) 
28(1) 
4(1) 

6(6) 
83(1) 
11(11) 

29(29) 
56(7) 
15(15) 



1(2) 

13(17) 

80(19) 

5(10) 

2(4) 

23(3) 

26(14) 
21(14) 

3(7) 
27(12) 

47(2) 

0(0) 
0(0) 

53(2) 

70(3) 

2(1) 
5(3) 

1(1) 
10(3) 
12(3) 

6(2) 
94(2) 



22(1) 

29(6) 
49(7) 

9(2) 
25(6) 
4(3) 
62(6) 



41(4) 

5(7) 

33(10) 

16(8) 
6(9) 

8(3) 
28(9) 
64(9) 

23(8) 
67(10) 
10(1) 



0-8 
5-15 
0-35 



40-50 
1-3 

50-60 
60-70 

1-5 
0-12 

0-2 
6-14 
5-20 



10-20 



20-30 
7-25 

30-60 



10-30 



15-30 



P33(1232) 



ttN 



100(1) 



100(0) 



98-100 



P33(1600) 



P33(1920) 



D33(1700) 



D35(1930) 



D35(2350) 



F35(1752) 



F35(1905) 



F37(1950) 



(7rA)p 
7rN*(1440) 

ttN 

{■kA)p 

7rN*(1440) 

ttN 

(7rA)p 

7rN*(1440) 

ttN 
(p3N)s 
{■kA)d 
{■kA)s 

ttN 
KA 
Flux 

ttN 
KA 
Flux 

ttN 

(P3N)P 

{nA)F 
{nA)p 

ttN 

(P3N)P 

(7rA)p 
(7rA)p 

ttN 

(p3N)p 

(7rA)p 
Flux 



0(1) 
0(1) 

28(5) 
59(10) 
13(4?1) 

5(61?) 
41(2.9) 
53(8.2) 

5(1.6) 

1(1) 

4(1) 

90(1.7) 

9(8) 
91(11) 



7(14) 
93(15) 



0(1) 
60(60) 
40(1) 

0(1) 

9(2.2) 
24(1) 
44(1) 
23(1) 

44(1) 

36(1) 
20(20) 



0(0) 
0(0) 

12(2) 
67(5) 
20(4) 

2(2) 
83(26) 
15(24) 

14(6) 
8(4) 
4(3) 
74(7) 

18(2) 

82(2) 

2(0) 

98(0) 

2(1) 
22(14) 
48(16) 
28(18) 

12(3) 
86(3) 
0(1) 
1(3) 

38(1) 
43(1) 
18(3) 



10-25 
40-70 
10-35 

5-20 



10-20 
5-20 
1-7 

25-50 

5-20 



5-15 
0-60 
0-25 



35-40 
0-10 
20-30 



Sii( 


1535) 


Sii( 


1650) 


Sii( 


2090) 


Pll 


1440) 


Pll 


1710) 


Pll 


2100) 


Pl3 


1720) 


Di3 


1520) 


Di3 


'1700) 


Di3 


'2080) 


Di5 


1675) 


Fl5 


1680) 


Fl7 


1990) 


Gl7 


;2190) 


S3l( 


1620) 


S3l( 


1900) 


P31 


1750) 


P31 


1910) 


P33 


1232) 


P33 


1600) 


P33 


1920) 


D33 


'1700) 


D35 


1930) 


D35 


;2350) 


F35 


1752) 


F35 


1905) 


F37 


1950) 



TABLE IX. Pole positions. The complex energy of the pole for each state is given along with the physical mass. 

Res Mass Pole Position 

(MeV)) (MeV) 

1545 1525 - 51 i 

1693 1663 - 120 i 

1822 1795 - 110 i 

1479 1383 - 158 i 

1699 1679 - 66 i 

2083 1810 - 311 i 

1716 1692 - 47 i 

1520 1504 - 56 i 

1729 1704 - 78 i 

2002 1824 - 307 i 

1687 1674 - 60 i 

1679 1667 - 61 i 

2311 2301 - 101 i 

2168 2107 - 190 i 

1633 1607 - 74 i 

1798 1795 - 29 i 

1721 1714 - 34 i 

1995 1880 - 248 i 

1234 1217 - 48 i 

1687 1599 - 156 i 

1889 1880 - 60 i 

1732 1726 - 59 i 

1932 1883 - 125 i 

2459 2427 - 229 i 

1724 1697 - 56 i 

1873 1793 - 151 i 

1936 1910 - 115 i 



FIG. 1. Schematic diagram for the Dyson equation iteration. 
FIG. 2. Schematic diagram for the self energy (S) iteration. Each resonance is allowed to couple to all open channels. 

FIG. 3. (ji{s) (channel propagator) distributions for the ttN and rjN channels in Sn partial wave. These functions are used 
in the two channel- two resonance test case and in the final fits. 

FIG. 4. One resonance-one channel test case T matrix (real and imaginary parts) using parameters appropriate for the 
P33(1232) resonance. Solid dots are placed at the energy of the resonance mass. 

FIG. 5. Self energy function, ^1535,1535 for the test case where only -kN and rjN channels are included. 

FIG. 6. T matrices for two channel-two resonance test case. The real and imaginary parts are shown for the nN — > nN and 
ttN -^ rjN reactions. Solid dots are placed at the values of the resonance masses obtained from these T matrices. 

FIG. 7. Best fit calculation for Sn partial wave at the i]N threshold. T matrices for nN — > ttN and nN -^ p\N are shown. 

FIG. 8. T-matrices for Sn partial wave for elastic scattering and final states of piA'^, crN, and 7rAf*(1440) calculated with 
the full model using final fit parameters. The dotted line is calculated with only resonance couplings turned on and the dashed 
line is calculated with only non-resonant couplings enabled. The full calculation, which cannot be a sum of the dotted and 
dashed lines because the resonant and non-resonant diagrams interfere, is shown as a solid line. 

FIG. 9. T matrices according to the best fit for the D15 partial wave. The VPI partial wave amplitudes are shown by data 
points with error bars. The KSU fit of Mauley and Saleski is shown as dashed lines and the fit of this work is shown as solid 
lines. 

FIG. 10. T matrices for D13 partial wave. The same labeling is used as in Fig. H. 

FIG. 11. T matrices for Sn partial wave. The same labeling is used as in Fig. H. 

FIG. 12. Total cross section for the n^p — > rjn reaction. The data was taken from sources in Ref. [27], leaving out the low 
energy Brown, et al. data. The solid line shows the total cross section calculated with the final fit result of this work using all 
partial waves. The contribution from the Sn partial wave would be very similar to the full result. 

FIG. 13. Pitt-ANL model calculation using final fit compared with the total cross section for the process n^p -^ Tv°n. The 
fit was to T-matrix amplitudes derived in part from the data shown here. Data shown is from the VPI SAID database. 

FIG. 14. Pitt-ANL model calculation using final fit parameters compared to the differential cross section for the process 
TT^p — ^ TT^p. Data from the VPI SAID database is shown. 

FIG. 15. Pitt-ANL model calculation using final fit parameters of the polarization, P, for the process 7r~p ^ n^p. Data 
from the VPI SAID database is shown. 
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